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Abstract 



We construct systematic expansions around four and two spatial dimensions for a Fermi gas near 
the unitarity limit. Near four spatial dimensions such a Fermi gas can be understood as a weakly- 
interacting system of fermionic and bosonic degrees of freedom. To the leading and next-to-leading 
orders in the expansion over e = 4 - d, with d being the dimensionality of space, we calculate the 
thermodynamic functions and the fermion quasiparticle spectrum as functions of the binding energy 
of the two-body state. We also show that the unitary Fermi gas near two spatial dimensions reduces 
to a weakly-interacting Fermi gas and calculate the thermodynamic functions and the fermion quasi- 
particle spectrum in the expansion over e = d-2. 

Then the phase structure of the polarized Fermi gas with equal and unequal masses in the unitary 
regime is studied using the e expansion. We find that at unitarity in the equal mass limit, there is 
a first-order phase transition from the unpolarized superfluid state to a fully polarized normal state. 
On the BEC side of the unitarity point, in a certain range of the two-body binding energy and the 
mass diff'erence, we find a gapless superfluid phase and a superfluid phase with spatially varying 
condensate. These phases occupy a region in the phase diagram between the gapped superfluid phase 
and the polarized normal phase. 

Thermodynamics of the unitary Fermi gas at finite temperature is also investigated from the per- 
spective of the expansion over e. We show that the thermodynamics is dominated by bosonic excita- 
tions in the low temperature region T T^. Analytic formulas for the thermodynamic functions as 
functions of the temperature are derived to the lowest order in e in this region. In the high temperature 
region where T ~ T^, bosonic and fermionic quasiparticles are excited and we determine the critical 
temperature and the thermodynamic functions around Tc to the leading and next-to-leading orders 
in 6 and e. 

Finally we discuss the matching of the two systematic expansions around four and two spatial 
dimensions in order to extract physical observables at J = 3. We find good agreement of the results 
with those from recent Monte Carlo simulations. 
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Chapter 1 



Introduction 



BCS-BEC crossover and the unitarity limit 

Interacting fermionic systems appear in various subfields of physics. The Bardeen-Cooper-SchriefFer 
(BCS) mechanism shows that if the interaction is attractive, the Fermi surface is unstable toward the 
formation of Cooper pairs and the ground state of the system universally exhibits the superfluidity or 
superconductivity mill- Such phenomena have been observed in the metallic superconductor |l3l, the 
superfluid ^He flU, the high-Tc superconductor JSl, and recently, in the ultracold atomic gases such 
as '*'^K |[6l or ^Li |[71 under optical traps. Possibilities of the superfluid nuclear matter |[8l|9l, the color 
superconductivity |[IOl[Ii|, and the neutrino superfluidity [12J are also discussed in literatures, which 
will be important to astrophysics such as the physics of neutron stars. Among others, experiments 
on the ultracold atomic gases have the remarkable feature where the strength of attraction between 
atoms is arbitrarily tunable through the magnetic-field induced Feshbach resonances [fT3l[T4| . Their 
interaction predominantly arises from binary 5-wave collisions whose strength can be characterized 
by the 5- wave scattering length a. Across the Feshbach resonance, a can in principal be tuned to take 
any value from -0 to -00 and from +00 to +0 lfT5l[T6l . Therefore, the experiments of the ultracold 
atomic gases provide an ideal field for studying interacting fermionic systems and testing various 
many-body techniques developed for related problems. 

With the use of the experimental technique of the magnetic-field tunable Feshbach resonance, the 
longstanding idea of the crossover from the BCS state to the Bose-Einstein condensation (BEC) ifTTl 
[HmHEDJ has been recently realized in laboratories. The schematic phase diagram of the BCS-BEC 
crossover problem is shown in Fig. 11.11 as a function of the dimensionless variable IKak^) with kp 
being the Fermi momentum (horizontal axis). Here the effective range of the interaction tq is assumed 
to be small tq^f 0. As long as the attractive interaction between fermions is weak (BCS regime 
where afcp < 0), the system exhibits the BCS ground state characterized by the condensation of the 
loosely bound Cooper pairs. On the other hand, if the attractive interaction is strong enough (BEC 
regime where < ak^), two fermions form a bound molecule and the ground state of the system is 
described by the BEC of the tightly bound molecules. These two apparently different situations are 
considered to be smoothly connected without the phase transition, which means the ground state in 
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^ I \ ak-p I =00 




Figure 1.1: The BCS-BEC crossover problem in the plane of the inverse scattering length l/ak^ and 
the spatial dimension d. There are four limits where the system becomes non-interacting; afcp ±0, 
d ^ 4, and d 2. 

the whole regime of the l/(akp) axis is the superfluidity/superconductivity. 

The both ends of the phase diagram correspond to weakly-interacting systems; a weakly-interacting 
Fermi gas in the BCS limit ak^ -0 and a weakly-interacting Rose gas in the EEC limit akp +0. 
In these two limits, the usual Green's function techniques for the many-body problem are applica- 
ble and we can understand the properties of the system in terms of the perturbative expansion over 
\akji\ «; 1 [ET]| . An exception is the superfluidity in the BCS limit where the paring gap A ~ e^/(2aiF) 
a non-perturbative effect and never appears to any finite orders in the perturbative akp expansion. On 
the other hand, there exists a strongly-interacting regime in the middle of the phase diagram where 
lakji] » 1. The point of the infinite scattering length ja^pl ^ 0° corresponds to the threshold where 
the zero energy bound molecule is formed and provides the conceptual boundary between the BCS 
regime and the BEC regime. The perturbative expansion over [ak^l obviously breaks down in this 
regime because the expansion parameter is no longer small. Such a regime is frequently referred to as 
the unitary regime because the infinite scattering length \akp\ ^ 00 is the "unitarity limit" where the 
5-wave cross section reaches the maximal value allowed by the unitarity constraint of the scattering 
matrix. 

The Fermi gas in the unitarity limit, which we refer to as the unitary Fermi gas, has attracted 
intense attention across many sub fields of physics. As we mentioned, the system can be exper- 
imentally realized in atomic traps using the Feshbach resonance and has been extensively stud- 
ied [|2ai23IMl25llMl27llMl29l[3a[3n[a[3S 

Furthermore, since the fermion density is the only dimensionful scale of the unitary Fermi gas, its 
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properties are expected to be universal, i.e., independent of details of the interparticle interaction. 
The unitary Fermi gas is thus an idealization of dilute nuclear matter f4l], where the neutron-neutron 
i'-wave scattering length a„„ =^ -18.5 fm is much larger than the effective range ro - 2.7 fm [48], and 
may be relevant to the physics of neutron stars. It has been also suggested that its understanding may 
be important for the high-Tc superconductivity H9l . The Fermi gas in the unitary regime provides 
a new interesting regime other than the conventional weak-coupling BCS and BEC regimes and is 
worthwhile to study experimentally and theoretically. 

Experimental achievements 

Before going to details of the theoretical treatment of this thesis, we briefly summarize related experi- 
mental achievements in the trapped fermionic atoms over the past decade. Typically 10^ ~ 10^ atoms 
are confined in the optical potential and can be cooled down to the temperature ~ 0. 1 using the 
laser cooling and evaporative cooling lISOl . where is the Fermi energy of order micro-Kelvin. The 
degenerate Fermi gas of atoms was first reported in [[STl, where the non-classical momentum distribu- 
tion and total energy of the gas were observed. The magnetic-field induced Feshbach resonance make 
it possible to realize the degenerate Fermi gas in the strongly-interacting regime where |a|^F ^ 1 • Ex- 
pansion dynamics of such a strongly-interacting degenerate Fermi gas was studied and an anisotropic 
expansion was observed when the gas is released from the optical trap E^ . The hydrodynamics of 
expanding strongly-interacting fermionic systems has also attracted special interest of nuclear physics 
community in connection with the strongly-interacting matter created in the experiment of Relativistic 
Heavy Ion Collisions JSl [53l [MlESl [51 [53 • 

The next goal after creating the degenerate Fermi gas of atoms in the strongly-interacting regime 
was to establish the superfluidity. First a long-lived molecular gas was created from the atomic Fermi 
gas using the Feshbach resonance |[25l l26l l27l l28l [32l and the condensation of the weakly bound 
molecules was observed on the BEC side a/tp > [|29l [301 1211 [331 [36l llQjl. The BEC of the bound 
molecules formed by fermionic atoms is in complete analogy with that of bosonic atoms achieved 
earlier using *^'Rb or ^^Na Ii58ll59l . The sudden onset of a peak in the momentum distribution of the 
molecules was reported indicating the phase transition to the BEC state. Soon after the realization of 
the molecular condensates on the BEC side, the condensation of fermionic atom pairs was observed 
on the BCS side akp < ^U^. Here a pairwise projection of the fermionic atoms onto molecules 
by a rapid sweep from negative a to positive a was used to measure the momentum distribution of 
the fermionic atom pairs. As an evidence of the superfluidity, collective excitations of the trapped 
Fermi gas were studied, where the measured frequencies and damping rates are plausibly explained 
assuming the superfluidity [|34l [3511381 . Another evidence of the superfluidity was an appearance of 
the pairing gap measured using the radio frequency spectroscopy [37] . The critical temperature as an 
onset of the superfluidity was presented from the measurement of the heat capacity [i39il . Finally the 
definitive evidence of the superfluidity was provided by the observation of vortex lattices in a rotating 
Fermi gas [|4T|. 
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The striking properties of these trapped atomic systems lay not only on the tunability of the in- 
teraction strength between atoms but also on their high designability of various interesting systems. 
For example, Fermi gases with a population imbalance between the two fermion components (finite 
polarization) have been already realized and studied fi421 H31 HH HH |46l . The phase separation in the 
trapping potential was observed where the superfluid core with equal densities is surrounded by the 
normal gas with unequal densities. Also the realization of Feshbach resonances is available between 
two different atomic species [|60ll6T]| . While such Feshbach resonances are currently achieved only 
between fermionic and bosonic atoms, it may be possible in future to realize the Feshbach resonances 
between two different species of fermionic atoms. It will bring us a new ground to investigate a 
Fermi gas with finite mass difference between different fermion species. Such asymmetric systems 
of fermions with density and mass imbalances will be also interesting as a prototype of high den- 
sity quark matter in the core of neutron stars, where the density and mass asymmetries exist among 
diff"erent quark flavors ill EH [Ml ISl El • 

Theoretical treatment and e expansion 

The breathtaking experimental progresses with high tunability and designability have stimulated var- 
ious theoretical studies in the new regime of the unitarity limit. The austere simplicity of the unitary 
Fermi gas, however, implies great difficulties for theoretical treatment, because there seems to be no 
parameter for a perturbation theory. The usual Green's function techniques for the many-body prob- 
lem is completely unreliable here since the expansion parameter akp becomes infinite in the unitarity 
limit. Although the mean-field type approximations (with or without fluctuations) are often adopted to 
give a qualitative picture of the BCS-BEC crossover problem, they are not controlled approximations 
in the strong-coupling unitary regime [[l7l[I3[IlEni[67l[Ml[63|7ni|7D|71|73|^ 
challenging problem for many-body theorists is to establish a systematic approach to investigate the 
Fermi gas in the unitary regime. 

At very high temperatures, a systematic expansion in terms of the fugacity z = e^^^ (virial expan- 
sion) is applicable where /u denotes the chemical potential and T denotes the temperature 11791 [SOllSTIl . 
However, the virial expansion is restricted to the normal state and can not describe the phase transition 
to the superfluid state. Considerable progress on the study of the unitary Fermi gas has been recently 
made by Monte Carlo simulations both at zero temperature |[H21[H3[Hll[H5l[86l[H21 and finite temper- 
ature [l88l [89l l90l l9n l92l . but these simulations also have various limitations. For example, Monte 
Carlo simulations can not treat systems with a population or mass imbalance between two fermion 
species due to fermion sign problem. Also they can not answer questions related to dynamics like the 
dynamical response functions and the kinetic coefficients. Analytic treatments of the problem at any 
temperature, if exist, will be extremely useful and give insights that are not obvious from numerics. 

Recently we have proposed a new analytic approach for the unitary Fermi gas based on the sys- 
tematic expansion in terms of the dimensionality of space Ii93ll94ll95l . utilizing the specialty of four 
or two spatial dimensions in the unitarity limit [96] . In this approach, one would extend the problem 
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to arbitrary spatial dimensions d with keeping the scattering length to be infinite la/rpl ^ oo (vertical 
axis in Fig. 11.11) . Then we find two non-interacting limits on the d axis, which are d = A and d = 2. 
Accordingly, slightly below four or slightly above two spatial dimensions, the unitary Fermi gas be- 
comes weakly-interacting, where a "perturbative expansion" is available. Actually the unitary Fermi 
gas near four spatial dimensions can be understood as a weakly-interacting system of fermionic and 
bosonic quasiparticles, while near two spatial dimensions it reduces to a weakly-interacting Fermi 
gas. A small parameter of the perturbative expansion there is e = 4 - J near four dimensions or 
e = d - 2 near two dimensions. The weakness of the interaction near four and two spatial dimen- 
sions can be also seen from the viewpoint of the renormalization group fixed point; the fixed point 
describing the unitarity limit approaches to the non-interacting Gaussian fixed point near d = 4 and 
d = 2 [|97ll98]| . After performing all calculations treating e or 6 as a small expansion parameter, results 
for the physical case of three spatial dimensions are obtained by extrapolating the series expansions 
to 6(6) = 1, or more appropriately, by matching the two series expansions. 

The 6 expansion around four spatial dimensions has been developed to calculate thermodynamic 
functions and the fermion quasiparticle spectrum at zero temperature near the unitarity limit to the 
leading and next-to-leading orders in e. The results were found to be quite consistent with those 
obtained by the Monte Carlo simulations and the experiments [|93l l94l . Then the 6 expansion was 
extended to investigate the thermodynamics of the unitary Fermi gas at finite temperature T, below 
and above the critical temperature T = [|95l . It should be noted that the critical dimension of a 
superfluid-normal phase transition is also four, which makes weak-coupling calculations reliable at 
any temperature for the small e. The e expansion was also successfully applied to study atom-dimer 
and dimer-dimer scatterings in vacuum and found to give results quite close to non-perturbative nu- 
merical solutions at J = 3 [|99l . Thus there are compelling reasons to hope that the limit J ^ 4 is not 
only theoretically interesting but also practically useful, despite the fact that the expansion parameter 
6 is one at J = 3. We recall that the e expansion has been extremely fruitful in the theory of the sec- 
ond order phase transition HIOOL In addition to the above-mentioned works, the phase structure of the 
polarized Fermi gas near the unitarity limit has been investigated based on the e expansion [|94lll01L 
The next-to-next-to-leading order correction to the thermodynamic functions at zero temperature was 
computed in HI 0211 and the BCS-BEC crossover was studied by the e expansion [I103L Very recently, 
another systematic approach to the Fermi gas in the unitary regime has been proposed; the two- 
component Fermi gas is generalized to a 2A^-component Fermi gas with Sp(2A'^) invariant interaction. 
Then the inverse of the component number 1 /N is used for the small expansion parameter [|98l I104L 

In this thesis, we give a comprehensive study of the unitary Fermi gas from the perspective of the e 
expansion. We start with the study of two-body scattering in vacuum in the unitarity limit for arbitrary 
spatial dimensions 2 < d < 4, which will clarify why and how systematic expansions around four and 
two spatial dimensions are possible for the unitary Fermi gas (Chap.|2l). In Chap.|3l we give a detailed 
account of the e expansion for the unitary Fermi gas and show the results on the thermodynamic 
functions and the fermion quasiparticle spectrum in the unitary regime to the leading and next-to- 
leading orders in e. Then the e expansion is applied to the unitary Fermi gas with unequal densities 
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of two components, which we call "polarized" Fermi gas. The phase structure of the polarized Fermi 
gas as a function of the two-body binding energy is investigated with equal fermion masses in Chap.lH 
and with unequal fermion masses in Chap [51 We also show in Chap. [6] that there exists a systematic 
expansion for the Fermi gas in the unitarity limit around two spatial dimensions. In Chap. Ul we 
make an exploratory discussion to connect the two systematic expansions around four and two spatial 
dimensions in order to extract results at d = 3. Chaps. [8] and |9] are devoted to the thermodynamics 
of the unitary Fermi gas at finite temperature. In Chap. [81 analytic formulas for the thermodynamic 
functions in the low temperature region T <^ Tc are shown as functions of the temperature to the 
lowest order in e. The behavior of the thermodynamic functions above to the leading and next-to- 
leading orders in e is discussed in Chap. [9] In particular, we put emphasis on the determination of the 
critical temperature and the thermodynamic functions at Tc by matching the e expansion with the 
expansion around two spatial dimensions. All the results thus obtained are compared to those from 
the recent Monte Carlo simulations. Finally, summary and concluding remarks are given in Chap. [TOl 
The chapters [2]-[2] in this thesis are based on the author's works in collaboration with D. T. Son |[93l 
[941I105II . while the chapters [8] and [9] are based on the author's paper [[95l . 



Chapter 2 



Two-body scattering in vacuum 



The special role of four and two spatial dimensions in the unitarity limit has been recognized by 
Nussinov and Nussinov [|96l . They noticed that at infinite scattering length, the wavefunction of two 
fermions with opposite spin behaves like R(r) = \/r^~^ at small r, where r is the separation between 
two fermions. Therefore, the normalization integral of the wavefunction takes the form of 



which has a singularity at rg ^ in high dimensions d > A. From this observation, it is concluded 
that the unitary Fermi gas at J ^ 4 become a non-interacting Bose gas. In particular, at fixed Fermi 
momentum the energy per particle goes to zero as J ^ 4 from below. On the other hand, in low 
dimensions d < 2, the attractive potential with any strength causes zero energy bound states, and 
hence, the threshold of the appearance of the first two-body bound state corresponds to the zero 
coupling. It follows that the Fermi gas in the unitary limit corresponds to a non-interacting Fermi 
gas at d ^ 2 and the energy per particle approaches to that of the free Fermi gas in this limit. This 
singular character of J = 2 was also recognized in the earlier work [I106II . 

In this Chapter, we shall give precise meaning for their intuitive arguments at J = 4 and 2 from the 
perspective of diagrammatic approach. These two dimensions will provide foundations to construct 
systematic expansions for the unitary Fermi gas. Here we concentrate on the two-body problem. 

2.1 Near four spatial dimensions 

The system under consideration is described by a Lagrangian density with a local four-Fermi interac- 
tion (here and below ^ = 1): 



where mo- (= m^, nii) is the fermion mass and cq is the bare attractive coupling between two fermions. 
Here we consider the fermions with equal masses = = m and generalize our discussion to the 
unequal mass case in Chap. [51 




(2.1) 
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Figure 2.1: Two-fermion scattering in vacuum in the unitarity limit. The T-matrix near four spatial 
dimensions is expressed by the propagation of boson with the small effective coupling g, while it 
reduces to a contact interaction with the small effective coupling near two spatial dimensions. 



The T-matrix of the two-body scattering is given by the geometric series of bubble diagrams 
depicted in Fig. 12.11 As a result of the summation, its inverse can be written as 

1 



1 r dkodk 1 



+ iS + kr, - sp+u + i6 



\ r dk 1 
7o~ } Jln^ 



(2.3) 



where Sp = p^/2m is the kinetic energy of non-relativistic particles. The integral in Eq. (12.31) is 
divergent and needs to be regularized. In this thesis we shall work in dimensional regularization: 
integrals are evaluated for those values of d where they converge [for the integral in Eq. (|2.3I) this 
corresponds to J < 2] and then analytically continued to other values of d. In this regularization 
scheme, the integral in Eq. (12.31) vanishes for po = p = (scatterings at threshold). Therefore the 
limit of infinite scattering length where r(0, 0) = oo corresponds to Cq = oo. 
The integration over k in Eq. (12.31) can be evaluated explicitly 

This expression has a pole at J = 4 originating from the ultraviolet divergence of the k integration. 
Substituting d = 4-6 and expanding in terms of e, the T-matrix near four spatial dimensions becomes 

iT{po,p) = + O(e^). (2.5) 

If we define 

' (2.6) 



and 



D(po,p) = {po-f + i6] . (2.7) 



2.2. Near two spatial dimensions 
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the r-matrix, to leading order in e, can be written as 

inpo,p)^(igfiDipo,p). (2.8) 

The function D(po,p) is the propagator of a particle with mass 2m. It is natural to interpret this 
particle a bound state of two fermions at threshold. We will refer to this particle simply as the boson. 
Eq. (12.81) states that the two-fermion scattering near four spatial dimensions can be thought of as a 
process that occurs through the propagation of an intermediate boson, as depicted in Fig. 12. 1[ The 
effective coupling of two fermions into the boson is given by g. An important point is that the effective 
coupling g ~ e^^^ is small near four dimensions. This fact indicates the possibility to construct a 
perturbative expansion for the unitary Fermi gas near four spatial dimensions in terms of the small 
parameter e. 



2.2 Near two spatial dimensions 

Similarly, the perturbative expansion around two spatial dimensions is possible. This is because the 
inverse of T-matrix in Eq. (12.41) has another pole at J = 2. Substituting d = 2 + e and expanding in 
terms of e, the T-matrix near two spatial dimensions becomes 

iT(po,p) = i—e + 0{f). (2.9) 
m 

If we define the effective coupling at 2 + 6 dimensions as 

f = -e, (2.10) 

m 

the T-matrix to the leading order in e can be written as 

iT{po,p)^ig\ (2.11) 

We see that the T-matrix near two spatial dimensions reduces to that of a contact interaction with the 
small effective coupling ~ e as depicted in Fig. 12. 1[ In this case, the boson propagator D{p) in 
Eq. (12.81) corresponds to just a constant -1. We note that the same effective couplings near four and 
two spatial dimensions in Eqs. (12.61) and (12.101) can be obtained by the fixed point of the renormaliza- 
tion group flow describing the unitarity limit Il97ll98ll . We defer our discussion of the expansion over 
6 = J - 2 to Chap. [6] and concentrate on the expansion over e = A - d. 



2.3 Binding energy of two-body state 

We shall be interested not only in the physics right at the unitarity point, but also in the vicinity of it. 
In other words, we shall assume that 1/co can be nonzero in dimensional regularization. The case of 
Co < corresponds to the BEC side of the unitarity point, and co > corresponds to the BCS side. 
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To facilitate a comparison with the physics in three dimensions, we shall find the relation between the 
coupling Co, assuming that we are on the BEC side, co < 0, and the binding energy of the two-body 
state at arbitrary spatial dimensions 2 < d < 4. 

The binding energy s^, defined to be positive > 0, is obtained as the location of the pole of the 
T-matrix at zero external momentum: T(-£\„ 0)"' = 0. From Eqs. (|2.3I) and (12.41) . we see cq is related 
with Sb via the following equation 



d 

d 



-=ni-^(^r=r. (2.12) 



Co \ IJXAtiI 

Near four spatial dimensions, the relationship between eb and Co to the next-to-leading order in 6 = 
4 - J is given by 

In three spatial dimensions, the relationship between Sb and the scattering length a = -mco/4n be- 
comes eb = (ma^y^ . If we recall that the Fermi energy is related to kp through sp = k^/(2m), one 
finds the following relationship holds at J = 3: 



(2.14) 



We will use these relations to compare results in the expansion over e = 4 - d with the physics at 
three spatial dimensions. 



Chapter 3 

Unitary Fermi gas around four spatial 
dimensions 



3.1 Lagrangian and Feynman rules 

Employing the idea of the previous Chapter, we now construct the e expansion for the Fermi gas near 
the unitarity limit. We start with the Lagrangian density in Eq. (|2.2I) with chemical potentials /i^ and 
fil introduced to two different components: 

^ " Zi ('^' ^ ^ ^ ^ ^° 'Aj'Al'Ai'AT- (3.1) 

D-=T,i ^ ' 



After making the Hubbard-Stratonovich transformation ni07[ I108II . the Lagrangian density in the 
Nambu-Gor'kov formalism can be written as 

X = (idt + +iuo-3+h]i'- —f(P + ^V+¥0 + ¥V_¥0*, (3.2) 

\ 2m / Co 

where *F = (i/^-i-, i/rj)^ is a two-component Nambu-Gor'kov field, and (Xi 2,3 and a+ = ^(iTi + 20-2) are 
the Pauli matrices jl 10 1 1 109]. We define the average chemical potential as yu = (//j^ + //J/2 and the 
chemical potential difference as H = Qxf - /zj,)/2. 

The ground state at finite density system (at least when // = 0) is a superfluid state where 
condenses: (0) = (po with 0o being chosen to be real. With that in mind we expand (f> around its 
vacuum expectation value 0o as 

(f> = (f>o + gif, g = . (3.3) 

Here we introduced the effective coupling g ~ e^^^ in Eq. (12.61) . The extra factor {nKpojlny^^ was 
chosen so that ip has the same dimension as a non-relativistic field §. Since the Lagrangian density in 



'The choice of the extra factor is arbitrary, if it has the correct dimension, and does not affect the final results because 
the difference can be absorbed into the redefinition of the fluctuation field if. The particular choice of g in Eq. ( 13.31 ) will 
simplify the form of loop integrals in the intermediate steps. 
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, 1^^- , ► ^ ® -------H 

iG igc^- -igh/cQ -i^Q 

-> «^ X >^ ---«--^--« <- (X) ^ ^ 

iD if^B -ig(j)o/co -2ifi 

propagators £i vertices £2 vertices 

Figure 3.1: Feynman rules from the Lagrangian density in Eq. (13.41) . The two vertices on the last 
column come from X2, while the rest from Xi. Solid (dotted) lines represent the fermion (boson) 
propagator iG (iD). 



Eq. (|3.2I) does not have the kinetic term for the boson field (p, we add and subtract its kinetic part by 
hand. In other words, we rewrite the Lagrangian density as a sum of three parts, X = Xo + Xi + X2, 
where 

X„ = (id, + H+ + (r+<^o + c^-0o) ^ + ^* (idr + ^) ^ - ^ , (3.4) 



2m / \ 4m / cq 

Xi = ^YV+¥^ + ^¥V_¥^* + //^VgY + (ifi - —] (p*(p - —if - —if* , (3.5) 



Co/ Co Co 

£2 = -(p*Udt + — \(fi -2jd(p*(p . (3.6) 

As we shall soon see, 0o coincides, to leading order in e, to the energy gap in the fermion spectrum. In 
the unitarity limit, we have g^/co = 0. When cq is finite and negative, -g^/co - Sb gives the binding 
energy of boson to the leading order in e from Eq. (12.131) . Throughout this thesis, we consider the 
vicinity of the unitarity point where g^/co can be treated as a small quantity §. 

The part Xo represents the Lagrangian density of non-interacting fermion quasiparticles and bosons 
with the mass 2m, whose kinetic terms are introduced by hand in Xo and taken out in X2- The propa- 
gators of fermion and boson are generated by Xo- The fermion propagator G is a 2 x 2 matrix. 



{po + Hf-El + i6 



'po + H + Sp -00 

-00 po + H - Spj 



(3.7) 



where Ep = ^Sp + (f>^ is the excitation energy of the fermion quasiparticle. We use the i6 prescription 
where 6 = 0^ for positive po while 6 = 0~ for negative po. The boson propagator D is given by the 
same form as in Eq. (|2.7I) . 

Dipo,p) = i^Po-Y + i^) ■ (3-8) 



The precise meaning of the "small quantity" here is that g /cq compared to is of the order e, g /cq ~ e<po. 



3.2. Power counting rule of e 
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0(6 1) 


A'^ (fermion density) 




(Fermi energy) 


0(1) 


00 (condensate), A (energy gap) 




P (pressure), E (energy density) 


0(6) 


yu, jUb (chemical potentials) 



Table 3.1: Summary of powers of e in various quantities at zero temperature when the condensate 0o 
is regarded as 0(1). 



The parts Xi and Li generate the vertices which are depicted in Fig. 13.11 The first two terms 
in Xi represent the fermion-boson interactions, whose coupling is proportional to g and small in the 
limit 6^0. This coupling originates from the two-body scattering in vacuum in the unitarity limit 
which was studied in the previous Chapter (see Fig. l2.1l) . The third and fourth terms are insertions of 
chemical potentials, and the boson chemical potential 

ix^ = lix-^-, (3.9) 

Co 

to the fermion and boson propagators, respectively. We treat these two terms as perturbations because 
yu will turn out to be small, yu/0o ~ 6, and we limit ourselves to the vicinity of the unitarity point where 
we can regard g^/co to be so small that jUb/^o ~ f • The last two terms in L\ give tadpoles to the ip and 
(p* fields, which are proportional to -igcpolco. The condition of cancellation of tadpole diagrams will 
determine the value of the condensate 0o- 

Finally, we treat the two terms in X2 as additional vertices for the boson propagator as depicted in 
the last column of Fig. lSTTl which are proportional to -/Ho and -li^ where 

no(;?o,p) = Po-y- (3.10) 

These two vertices play roles of counter terms so as to avoid double counting of certain types of 
diagrams which are already taken into Xo and Xi . In the latter part of this Chapter, we will consider 
the unpolarized case with // = at zero temperature. 

3.2 Power counting rule of e 

Near unitarity, we can use the same power counting rule of e developed in the unitarity limit W^ . Let 
us first consider Feynman diagrams constructed only from Xo and Xi, without the vertices from X2- 
We make a prior assumption ///0o ~ f ^ which will be checked later, and consider 0o to be 0(1). Each 
pair of boson-fermion vertices is proportional to ~ 6 and hence brings a factor of e. Also, each 
insertion of yu ~ 6 or //b = 2/i - g^lc^ ~ e brings another factor of e. Therefore the naive power of e 
for a given diagram is Ngll + N^^, where Ng is the number of couplings g from Jl\. N^ = N^^ + Nf^^ is 
the sum of the number of fi insertions to the fermion line and /ub insertions to the boson line. 
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However, this naive counting does not take into account the fact that there might be inverse powers 
of e coming from integrals which have ultraviolet divergences at d = 4. Each loop integral in the 
ultraviolet region behaves as 



J dpodp ~ J 



dpsp-'p^ (3.11) 



while each fermion or boson propagator behaves, at worst, as G{p) ^ p ^ or D(p) ~ p ^. Therefore, 
a given diagram may diverges as ~ p^ with D being the superficial degree of divergence given by 

D = 6L-2Pp-2P^. (3.12) 

Here L is the number of loop integrals and /'f(b) is the number of fermion (boson) propagators. From 
the similar argument to the relativistic field theories Ulllil . one can derive the following relations 
including the number of external fermion (boson) lines £'f(b): 



(3.13) 



L = {Pf- N,,) + (Pb - N^,) -Ng + 1, 

N, = (Pf - A^^p) + Y = 2(P^ - N,J + £b. 

With the use of these relations, the superficial degree of divergence D is written as 

D = 6- 2(£f + Eb+ N^), (3. 14) 

which shows that the inverse powers of e appear only in diagrams with no more than three external 
lines and chemical potential insertions. This is similar to the situation in quantum electrodynamics 
where infinities occur only in electron and photon self-energies and the electron-photon triple vertex. 

However, this estimation of D is actually an over-estimate: for many diagrams the real degree 
of divergence is smaller than given in Eq. (|3.14l) . To see that we split G(p) into the retarded and 
advanced parts, G{p) = G^{p) + G^{p), where G^ (G^) has poles only in the lower (upper) half of the 
complex Po plane: 

Pq -|- £p Ep -\- Ep Ep Sp 

Gu(P) = —J —o = W7^^ ^ . -ex + 



- £ 2 + 2Ep (po -Ep + id) 2Ep (po + Ep- id) ' 

Q ( \ — Po ~ _ Ep ~ Ep + Sp 

''^^^ ~ p^ -EI + id ~ 2Ep (po -Ep + id) 2Ep (po + Ep- id) ' 
Gnip) = G2i(p) = —5 — ; = X77-^ ^ . -ex + 



p^ -El + id 2Ep (po -Ep + id) 2Ep (po + Ep- id) ' 

From these expressions, it is easy to see that the ultraviolet behaviors of different components of the 
propagators are different: 

Gfi(p) ~ GUp) - D\p) - p-\ (3.16) 

G\^\p) - Gl[\p) - p-\ GUp) - G^,ip) - p-'. 



3.2. Power counting rule of e 
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Note that the boson propagator D{p) has a pole on the lower half plane of po, we have only the 
retarded Green's function for boson, D^(p) = 0. 

From these analytic properties of the propagators in the ultraviolet region and the vertex structures 
in Xi as well as the relation of Eq. (|3.14l) . one can show that there are only four skeleton diagrams 
which have the l/e singularity near four dimensions. They are one-loop diagrams of the boson self- 
energy rFigs. l3.2r a) and l3.2r c)1. the (f tadpole rFig. l3.2r e)1. and the vacuum diagram rFig. l3.2r g)1. We 
shall examine these apparent four exceptions of naive power counting rule of e one by one. 



The first diagram, Fig. I3.2r a'). is the one-loop diagram of the boson self-energy. The frequency 
integral can be done explicitly, yielding 



2 r dk 



1 



X 



(3.17) 



£fc^^ + Ek+R -po- i5 Ek__R + £fc+^ +po- id 

The integral over k is ultraviolet divergent at J = 4 and has a pole at 6 = 0. Thus it is 0(1) by itself 
instead of 0(6) according to the naive counting. The residue at the pole is 

(3.18) 



na(p) = -g- 



-(po-y) + 0(6), 



which is cancelled out exactly by adding the vertex IIo in £.2- Therefore the diagram of the type in 
Figs. I3.2r a). when combined with the vertex from X2 in Fig. I3.2r b). conforms to the naive e power 
counting, i.e., is 0(e). 

Similarly, the diagram in Fig. I3.2r c) representing the boson self-energy with one p insertion is 
given by 



in,(,p) = 2pg- 



r dk 
J (2^ 



Gn{k + ^ 



G22U 



(3.19) 



which also contains a 1/6 singularity, and is 0(e) instead of the naive O(e^). The leading part of this 
diagram is 

dk 



U,ip) = -2pg- 



r dk 

J w 



(3.20) 



= -2p + 0{e\ 

and is cancelled out by the second vertex from X2- Then the sum of Figs. I3.2r c) and l3.2r d) is again 
0{6^), consistent with the naive power counting. 

The tadpole diagram with one p insertion in Fig. I3.2r e) is 

dk 



C d 
J (2;r 

r dk 

= igl^(f>o -pr- 
J (2;r 



dk su 
YlEl 



{Gn{k)-G22ik)] G2i{k) 



+ Oie"\ 



(3.21) 
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= 0(e3/2) 



Figure 3.2: Four apparent exceptions of naive power counting rule of e, (a, c, e, g). The boson self- 
energy diagram, (a) or (c), is combined with the vertex from X2, (b) or (d), to restore the naive e 
counting. The tadpole diagram in (e) cancels with other tadpole diagrams at the minimum of the 
effective potential. The vacuum diagram (g) is the only exception of the naive power counting rule of 
6, which is 0(1) instead of 0(6). 



which is 0{€^^^) instead of the naive 0(e^^^). This tadpole diagram should be cancelled by other 
tadpole diagrams of order 0(€^^^) in Fig. l3.2r f). 



—G2iik)-i— 

J ° (3.22) 

ig (m(pQ\^ ,g(pQ 



r d 



\ In I Cn 



2 \ In I Co 

The condition of cancellation, Sg + Sf = 0, gives the gap equation to determine (poip) to the leading 
order in e. The solution to the gap equation is 

0o = ^--(-f + O(.). (3.23) 



Co \ m / 

When cq < 0, this solution can be written in terms of the binding energy £b as cpQ = {2fi + Sb)/€. Now 
the previously made assumption fi/(f>o = 0(e) is checked. The condition of cancellation of tadpole 
diagrams is automatically satisfied by the minimization of the effective potential as we will see later. 

Finally, the one-loop vacuum diagram with one fi insertion in Fig. I3.2r g) also contains the 1/e 
singularity as 

dk Sk 



. ^ (3.24) 



3.3. Effective potential to leading and next-to-leading orders 
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0(1) 



0(1) 



0(e) 



Figure 3.3: Vacuum diagrams contributing to the effective potential up to the next-to-leading order in 
€. The second diagram is 0(1) instead of naive 0(e) because of the l/e singularity. 

The leading part of this diagram is 0(1) instead of naive 0(e) and can not be cancelled by any other 
diagrams. Therefore, Fig. l3.2r g) ~ 0(1) is the only exception of our naive power counting rule of e. 

Thus, we can now develop a diagrammatic technique for the Fermi gas near the unitarity limit 
where g^/\co\ ~ /i in terms of the systematic expansion of e = 4 - The power counting rule of e is 
summarized as follow: 

1. We consider fJ./(po ~ e and regard 0o as 0(1). 

2. For any Green's function, we write down all Feynman diagrams according to the Feynman rules 
in Fig. 13. II using the propagators from Xo and the vertices from Xi. 

3. If there is any subdiagram of the type in Fig. I3.2r a) or Fig. I3.2r c). we add a diagram where the 
subdiagram is replaced by a vertex from X2, Fig. I3.2r b) or Fig. I3.2r d). 

4. The power of e for the given Feynman diagram will be 0(6^^'^^^^''), where Ng is the number of 
couplings g and A''^ is the number of chemical potential insertions. 

5. The only exception is the one-loop vacuum diagram with one // insertion in Fig. I3.2r g). which 
is 0(1) instead of the naive 0(e). 

We note that the sum of all tadpole diagrams in Figs. I3.2r e) and l3.2r f) vanishes with the solution of 
the gap equation ^q. 

3.3 Effective potential to leading and next-to-leading orders 

We now perform calculations to leading and next-to-leading orders in e, employing the Feynman rules 
and the e power counting that have just been developed. To find the dependence of (()q on /u, we use 
the eff"ective potential method nill| in which this dependence follows from the the minimization of 
the eff"ective potential yeff(0o)- The eff"ective potential at the tree level is given by 




(3.25) 
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Up to the next-to-leading order, the effective potential receives contributions from three vacuum dia- 
grams drawn in Fig. 13.31 fermion loops with and without a pi insertion and a fermion loop with one 
boson exchange. 

The contributions from two one-loop diagrams are 0(1) and given by 



= ij j;^Tr[\nG-\p)+ncT,Gip)] 
dp 



(3.26) 



By changing the integration variable to z = {Spjcpof and with the use of the formula 

z"-^ Y{a)Y{p - a) 



r^OO 

Jo \z^lf 



(3.27) 



we can perform the integration over p that results in 

\ 2n ) 



Vi(0o) = - 



2r(f) 



d 

' T(-)r(-- - -) . , , 

00 - 



(3.28) 



r(-i) r(i) 

Substituting d = 4 - e and expanding in terms of e up to 0(e), we obtain the effective potential from 
the one-loop diagrams as 



Vi(0o) = Y 



7-3(y + ln2) 

1 + Z 6 



m(f)o 
In 



dl2 



l-2(y-ln2) 
1 + : e 



In 



dl2 



+ 0{e\ (3.29) 



where y ~ 0.57722 is the Euler-Mascheroni constant. 

The contribution of the two-loop diagram is 0{e) and given by 



V2(0o) = J Tr {G{p)(T,G{q)cT:\ D{p - q) 

2 r dpdq 
J (2;r)2^+2 



Gn{p)G22{q)D{p - q). 

Performing the integrations over pq and q^, we obtain 

Cdpdq {Ep- Ep){Eq- Sq) 



V2(<^0) 



= _£_ r ^ 

4 J {2n) 



(3.30) 



(3.31) 



T)^'^ EpEq{Ep Eq -\- Sp-q/2,) 

Since this is a next-to-leading order diagram and the integral converges at J = 4, we can evaluate it at 
d = 4. Changing the integration variables to x = £p/(f>o, y = ^q/(po, and cos 6 = p ■ q, the integral can 
be expressed by 



V2((po) = -e 



m<p(i 
In 



)dl2 poo poo pn 

— dx \ dy \ dOxysin^e 
^ Jo Jo Jo 



Vm-x\U{y)-y] 
f(x)f(y) [g(x,y) - V^cos6'] 



(3.32) 



with f{x) = Vx^TT and g{x,y) = f{x) + f(y) + ^(x + y). The integration over can be performed 
analytically to lead to 



V2m = -Cel^- 



2n 



00, 



(3.33) 



3.4. Thermodynamic quantities near unitarity 
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where the constant C is given by a two-dimensional integral 

[f(x)-x][f(y)-y] 



C = I dx I dy ■ 

Jo Jo mf(y) 

The numerical integrations over x and y result in 

C« 0.14424. 



[g(x,y) - ^Jg\x,y)-xy] . 



(3.34) 



(3.35) 



Now, gathering up Eqs. (13.251) . (13.291) . and (|3.32l) . we obtain the effective potential up to the next- 
to-leading order in e. 



Kff(0o) = Vo(0o) + Vi(0o) + 



(3.36) 



Co 



-I- 



|L 7-3(r.ln2)^_3 l_W^^l-2(r-ln2)^ 



2;r / 



The condition of the minimization of the effective potential in terms of 0o gives the gap equation; 
dVes/d(f>o = 0. The solution of the gap equation satisfies 



^„ = ai;[,H3c-i+i„2).]-?*£(^r 

€ Co \m I 



y -h In 2 
1 + 3C- 1 + ^- \e 



(3.37) 



Using the relation of cq with the binding energy of boson in Eq. (12.121) . we can rewrite the solution 
of the gap equation in terms of e\, up to the next-to-leading order in e as 



00 = ^ [1 + (3C - 1 + ln2) 6] + — 
6 e 



l+|3C-^ + ln2-^ln^|6 

2 2 00 



(3.38) 



We note that the leading term in Eq. (13.381) could be reproduced using the mean field approximation, 
but the 0{e) corrections are not. The 0(e) correction proportional to C is the result of the summation 
of fluctuations around the classical solution and is beyond the mean field approximation. 



3.4 Thermodynamic quantities near unitarity 

The value of the effective potential Vgff at its minimum determines the pressure P = -Veff(0o) at a 
given chemical potential pi and a given binding energy of boson s^. Substituting the dependence of 0o 
on ju and s\j in Eq. (|3.38l) . we obtain the pressure as 



6 



,17 7 + ln2 



3£b 

400 



~2n 



(3.39) 



The fermion number density A^^ is determined by differentiating the pressure in terms of /i as 



Ofi 6 



, 1 -27 + 21n2 
1 + —. e 



m(po 



'ill 



(3.40) 
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Then the Fermi energy from the thermodynamic of free gas in d spatial dimensions is given by 



2n 

£f = 

m 



ir\- + i\N 



2/d 



00 



1 -ln2 
1 : e 



(3.41) 



Note that non-trivial dependences on e like ^fe and In 6 appear by taking N ^ to the power of 2/d. 

From Eqs. (13.381) and (|3.41l) . we can determine the ratio of the chemical potential and the Fermi 
energy {x/sp near the unitary limit as 



.3/2 



■expl 



tine 
8-26 



5 



e eXne e 
1 + ;t + — ; ;t In — 

2 4 2 ep 



(3.42) 



The logarithmic terms in the second line originates by introducing 0o = f'^^Sp to the Insb/^o term 
in Eq. (13.381) . The first term in Eq. (13.421) gives the universal parameter of the unitary Fermi gas 

^ = y-ls^L^Q as 



.3/2 



1 + 



6 In 6 



8 



-|3C--(l-ln2) 



1 



3/2 



1 



(3.43) 



+ — e^/2ine_o.0246e^/2 + --- . 
16 



Here we have substituted the numerical value for C « 0.14424 in Eq. (13.351) . The smallness of the 
coefficient in front of e'^'^ is a result of the cancellation between the two-loop correction and the 
subleading terms from the expansion of the one-loop diagrams around J = 4. The 0(6^^^) correction 
to the universal parameter ^ was recently computed to find the large correction 0.480 e'^^ 01021 . This 
may be related with the asymptotic nature of the e expansion and some sort of resummation will be 
necessary to go beyond the next-to-leading order level. 

Using Eqs. (13.391) . (13.401) . and (|3.41l) . the pressure near the unitarity limit normalized by e-^N is 
given by 

(3.44) 



epA^ 



d^l 8ep^' 



Then the energy density E = [xN - P can be calculated from Eqs. (13.421) and (13.441) as 

Ed' 



d + 2 2ep 



6 6ln6 6, fib 

1 + T + — ; ;t In — 

4 4 2 ep 



(3.45) 



The pressure and energy density in the unitarity limit are obtained from ^ via the universal relations 
depending only on the dimensionality of space. Partial resummation of logarithmic terms of the 
binding energy ln(eb/£p) change the exponent of Sb/^F to (sb/eFY'^^'^ II103L 



3.5 Quasiparticle spectrum 

The e expansion we have developed is also useful for the calculations of physical observables other 
than the thermodynamic quantities. Here we shall look at the dispersion relation of fermion quasi- 
particles. To the leading order in e, the dispersion relation is given by a>p(p) = Ep = Js^ + (p^. 



3.5. Quasiparticle spectrum 
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Figure 3.4: Corrections to the fermion self-energy of order 0(e); a ju insertion and one-loop diagrams. 



which has a minimum at zero momentum p = with the energy gap equal to A = 0o = (2// + Sb)/e. 
The next-to-leading order corrections to the dispersion relation come from three sources: from the 
correction of in Eq. (13.381) . from the fi insertion to the fermion propagator, and from the one-loop 
self-energy diagrams, -iL{p), depicted in Fig. 13.41 

Using the Feynman rules, the one-loop diagram of the fermion self-energy in Fig. l3.4l is given by 



- n:{p) = / J [cT,G{k)cT_D{p -k) + cr_G(k)cr^Dik - p)] . (3.46) 
There are corrections only to the diagonal elements of the self-energy and each element is evaluated 



as 



^u(.p) 



= . 2 r dk 

^ J (2;r)^+i 



G22(k)D{p - k) 



_ r _dk Ek - 

2 J {2ny' Ek(Ek + Sk-p 



£fc 

12- pq- isy 



(3.47) 



and 



^22(P) = ig- j Gi,ik)D{k - p) 
~ 2 J (2n) 



(3.48) 



r)^ Ek(Ek + Sk-p/2 + po- i5) ' 



The dispersion relation of the fermion quasiparticle covip) is obtained as a pole of the fermion 
propagator det[G~H^^ P) + /^o-3 - ^(^^ P)\ = 0> which reduces to the following equation: 



(t>Q OJ + Sp-jU- 2:22(<^», p) 



0. 



(3.49) 



To find the 0(e) correction to the dispersion relation, we only have to evaluate the self-energy S(a;, p) 
with a> given by the leading order solution to = Ep. Denoting 'Zn(Ep,p) and 'Z22(Ep,p) simply by 
Sii and S22 and solving Eq. (13.491 ) in terms of to, we obtain the dispersion relation of the fermion 
quasiparticle as 

2ll - 222 -2fi . 2^ 



2i 1 -l- £22 
a>F{p) = Ep + + 



2Er 



-Sp + Oie'). 



(3.50) 
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cjf(p) 




A ~ 0(1) 

I' ^ 

>- — 

>-2 2m 

£o ~ 0(e) 



Figure 3.5: Illustration of the dispersion relation of the fermion quasiparticle in the unitarity limit 
from the e expansion, copip) = sjisp - sof- + in Eq. (13.531) is plotted as a function of Sp = 12m. 
Values on both axes are extrapolated results to 6 = 1 in units of the chemical potential yu. 



Since the minimum of the dispersion relation will appear at small momentum Bp ~ p, we can expand 
SC^p, p) around zero momentum p = as 2(£'p, p) = iPicpo, 0) + 'L'{(f)Q, 0) Sp/cpQ to find the energy 
gap of the fermion quasiparticles. Performing the integration over fc in 2 analytically, we have 

Sii(£p,p) = 6(2-81n3 + 81n2)0o + 6|-^ + 81n3-81n2j£p, (3.51) 

and 

^i2{Ep,p) = 6(-2-81n3 + 161n2)0o + 6|-^ -81n3 + 161n2|ep. (3.52) 



Introducing these expressions into Eq. (|3.50l) . we find the fermion dispersion relation around its min- 
imum has the following form 

co,(p) - A + - J(ep-eo)2 + A2. (3.53) 

Here A is the energy gap of the fermion quasiparticle, which is given by 

^ 2 (3.54) 

= [l -(8 In 3- 12 In 2) 6 + 0(6^)] (^0- 

The minimum of the dispersion curve is located at a nonzero value of momentum, |p| = (2meo)'^^, 
where 

, ^22 ^11 ^11 ^ ^22 

So=P + 

^ (3.55) 

= p^f + 0(eh 



3.6. Location of the splitting point 
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Then, introducing the solution of the gap equation (13.381) . the energy gap A as a function of the 
chemical potential and the binding energy is given by 



A = — [1 +(3C - 1 -81n3 + 131n2)6] + — 

6 6 



1 +I3C - ^ -81n3 + 131n2- ^In — |6 
2 2 A I 



= ^ (1 _ 0.345 .) + (l + 0.155 . - ln( 1 + ^ I , (3.56) 



e e \ I j I \ Eb 

while eo is given by 

£0 = 2// + ^. (3.57) 

Note the difference with the mean field approximation, in which eo = When sq is positive, the 
fermion dispersion curve has its minimum at nonzero value of momentum as in the BCS limit, while 
the minimum is located at zero momentum when eo is negative as in the BEC limit. We find the former 
(eo = 2/i > 0) is the case in the unitarity limit. The dispersion curve of the fermion quasiparticle ol)y:{p) 
in the unitarity limit eb = is illustrated Fig 13.51 In particular, the difference between the fermion 
quasiparticle energy at zero momentum and at its minimum is given by Sq/2^o = efx ~ 0(6^). 



3.6 Location of the splitting point 

Let us consider the situation where the binding energy eb is increasing from zero while the number 
density is kept fixed. Then A > is held fixed but // is changing. Since to the leading order in e, 
the chemical potential as a function of the energy gap is given by Iji = eA- eb, the location of the 
minimum of the dispersion curve can be written as 

£o = eA-^. (3.58) 

Therefore, eo decreases as eb increases. When the binding energy reaches eb = 26A, the minimum 
of the dispersion curve sits exactly at zero momentum p = 0. This point is referred to as a splitting 
point II112L We find the splitting point is located at the BEC side of the unitarity limit where the 
binding energy is positive eb = 2eA > and the chemical potential is negative 2p = -eA. This 
splitting point will play an important role to determine the phase structure of the polarized Fermi gas 
in the unitary regime as we will study in Chap.lH 



3.7 Momentum distribution function 

Other interesting observables which can be measured in experiments are a momentum distribution 
function of fermion quasiparticles 111131 1114[ 11151 1116II and a condensate fraction in the fermion 
density 1161171 11161 II 17[[TT^ . The momentum distribution functions, npf and npi, can be computed 
from the fermion propagator ^(p) through 

= J ^e'P'''gn{po,p) and n^^ = - J ^ e^'^'^''^ g22(Po,p), (3.59) 
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where the fermion propagator up to the order e is given by 

0(p) = G{p) + G{p) [I.{p) - pa^] G(p) + O(e^). (3.60) 

With the use of the definition of the bare fermion propagator G(p) in Eq. (13.71) and the one-loop 
fermion self-energy in Eq. (13.461) . the diagonal elements of the fermion propagator becomes 

Quip) = Gn(p) + Gnip) [I-n(p) - p] Gn(p) + Gn{p) [^22{p) + p] Gnip) 



Po + Sp 



+ 



p^-E^ + i6 {p^-E^ + id 



Po + Sp 



[Zuip)-p] + 



00 



p^-E^ + i6) 



^2 (3.61) 
[i:22ip)+p] 



and Q22{po, P) = -Qu(-po,p)- Since we are considering a symmetric case with p-^ = pi and = mi, 
we can confirm the relationship np = np^ = Upi. 

From the first term in Eq. (13.611) . the 0(1) part of the momentum distribution function rip^ is easily 
obtained as 

< = r ^ e'-^^'-p^^^ = (3.62) 



2Er 



'0 ' ^^P 

The last two terms in Eq. (13.611) contribute to the momentum distribution function to the order of e. 
Since 2ii(/>o,p) in Eq. (|3.47l) is analytic at Im[po] > 0, the po integral of the second term can be 
evaluated as 



dpo 
2m 



Po + Sp 



d 

[T.n{p)-p]= — 
opo 



Pa + Sp 
Po - Ep 



[^n{p)-p] 



\I.n{-Ep,p)- p] + 



po=-Ep 

(.Ep ~ Sp) 
~^Ef~ 



(3.63) 



2'ii(-£p,p). 



The prime in l,[^{pQ,p) represents the derivative with respect to po. Similarly, using the relationship 
^22(^0^ P) = -2ii(-Po, P), the Po integral of the third term in Eq. (13.611) can be evaluated as 

\2 ^ , / , \2 



r dpo 

J 2m 



00 



- E^ + iS 
V "0 p 



[l22(;.)+A^] = -/f^ 



01 







d 
dpo 



^pl-El + i6 

00 



[^n{p)-tA 



Po - Ep 



\Ln{-Ep,p) - p] 



[^u{p)-iA 
00 



(3.64) 



P0 = 



AEl 



I.[ii-Ep,p). 



From these two contributions, we obtain the next-to-leading order correction to the momentum dis- 
tribution function as 

1.2 



2EI 



2E 



\i:„(-^„p) - f^%^2'„(-£„p). 



p 



2El 



(3.65) 



Substituting the solution of the gap equation in the unitarity limit p = e(f)o/2+0{e^) and the expression 



for 'Lnip) in Eq. (13.471) . the Oie) part of the momentum distribution function rip^ is given by 



AEl 4 J {27if EpEuiEp + 



Sk) 



Ek + Sk-p/2) 



+ 



0i 

£2 ■ Ep(Ek + Sk-p/2 + Ep) 



(3.66) 



3.8. Condensate fraction 
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Figure 3.6: Momentum distribution function Hp of the fermion quasiparticles in the unitarity limit as 
a function of the momentum \p\/ yjlmcpQ to the leading and next- to-leading orders in e. The 0(1) part 
n^p (red curve), the 0{e) part Up^ (green curve), and their sum Hp = n^p + r{p (blue curve) at 6 ^ 1 
are respectively shown. 



The momentum distribution function of the fermion quasiparticles Up = rip'' + to the leading 
and next-to-leading orders in e is plotted as a function of the momentum \p\l -sjlm^o in Fig. l3.6[ Since 
the condensate 0o is related to the Fermi energy = k^/2m by Eq. (13.411) . 

00 /. l-ln2 \ 
^F = ^ 1 z e\, (3.67) 



the Fermi momentum in the horizontal axis of Fig. 13. 61 corresponds to 



_^-^-i/^(l_LJ!l^^|« 0.962 (6^1). (3.68) 



We can find overall agreement of the momentum distribution function computed by the e expansion 
up to the next-to-leading order in e with that from the Monte Carlo simulation II116II . 



3.8 Condensate fraction 



The fermion density in the condensation A^o is related with the off-diagonal element of the fermion 
propagator ^i2(;?) through lfTT6l [TTTl [TT8l [TT9ll 



(3.69) 
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Since Qn{p) up to the order e is given from Eq. (|3.60l) as 

Qn{p) = Guip) + Gnip) [^nip) - p] Gi2(p) + Gnip) [^iiip) + p] Gxiip) 



00 



00 (Pq + gp) 

pl-El + i6 (^p2_E2 + isy 

the po integration of each terms results in 

^dpo 



[^nip)-t^] 



00 



and 



r dpo 00 {Po + ^p) rv ^ N 1 r 



2;Ti p^- + id 

dpo 00 (Po - Sp) 



00 (Po - gp) 
(i^o' - + '■^)' 

00 



[S22(;?)+//], 



(3.70) 



(3.71) 



2ni 
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-00 

opo 



(pI - £ 2 + /5) 



[Z22(rt+y"] 



(po - Ep) 



[^n{p)-p] 



(3.72) 



= -00-^ [2n(-£p,p) + 0o%TT^2;i(-£p,p). 



Here we used the relationship 222(i?o,p) = -^n{-pii,p) and the fact that Sn(po,p) in Eq. (|3.47l) is 
analytic at Im[/)o] > 0. 

Therefore, the fermion density in the condensation A^o up to the order e is given by 



dp 

Y 



1 



2£2 £4^ £4 



Ep 



(3.73) 



p p p 

The first term provides the leading contribution to Nq, while the other three terms are 0{e) corrections. 
Substituting the solution of the gap equation in the unitarity limit p = 60o/2 + O(e^) and performing 
the integration over p, we obtain A^o as 



A^o = 



r(i-f)r(|) . 



+ -e + 0.127741 6 + 0(6^) 
8 



m0o 
2n 



(3.74) 



4r(2-f) 

Note that because the p integration of the leading term is logarithmically divergent at J = 4, the 
leading contribution to A^o becomes 0(1/6) at J = 4-6. The fermion density in the condensation 
should be compared to the total fermion density A^ to the same order in e II102II : 



m0o 
2n 



r(i-f)r(i + f) 1 

^ = _vi — l^L_v — ^ + + 0.258352 e + 0{e^) 
2V^r(2-f) 2 

Taking the ratio of A^o to A^, we find the condensate fraction Nq/N to be 

^ = 1 - 0.0966 6 - 0.2423 + 0{e^). 

The condensate fraction can take a value from zero (BCS limit) to one (BEC limit). The lowest term 
in the e expansion is one as in the BEC limit and thus all fermion are in the condensation at four 
spatial dimensions. Below d = 4, the condensate fraction decreases and the naive extrapolation to 
6 = 1 gives M,/A^ ~ 0.661. 



(3.75) 



(3.76) 



3.9. Extrapolation to e = I 
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A/lu 


No/N 


mean-field approximation 


0.5906 


1.1622 


0.6994 [|118|| 


6 expansion (NLO) 


0.475 


1.31 


0.661 


1/A^ expansion (NLO) 11041 


0.279 


1.50 




self-consistent approach [TTSlI 


0.360 






Monte Carlo simulation 


0.42 If86l 


1.2 il 


0.58 IIII6II 



Table 3.2: Comparison of the results by the e expansion with other analytic and numerical calculations 
in the unitarity limit. 



3.9 Extrapolation to 6 = 1 

Finally, we discuss the extrapolation of the series expansion to the physical case at three spatial 
dimensions. Although the formalism is based on the smallness of e, we see that the corrections are 
reasonably small even at 6 = 1 . If we naively use only the leading and next-to-leading order results 
for ^ in Eq. (l3.43L A in Eq. (l3.56L so in Eq. (l3.57L and Nq/N in Eq. (13.761) in the unitarity limit, their 
extrapolations to 6 = 1 give for three spatial dimensions 

^-0.475, -~1.31, — «2, —« 0.661. (3.77) 
jj. lu N 

They are reasonably close to the results of recent Monte Carlo simulations, which yield ^ « 0.42, 
A/// « 1.2, so/fi « 1.9 [[86l . and No/N ~ 0.58 [I116L They are also consistent with recent experimen- 
tal measurements of ^, where ^ = 0.51 + 0.04 [|39l and ^ = 0.46 + 0.05 [43 J. These agreements can 
be taken as a strong indication that the e expansion is useful even at 6 = 1 . In Table 13. 2[ we show 
results on ^, A/ep, and Nq/N from other analytic and numerical calculations. The self-consistent 
approach fTSl quoted in Table 13.21 and Table 19.11 is based on the Luttinger-Ward and DeDominicis- 
Martin formalism [|120[ I121| where the potential functional is self-consistently approximated by lad- 
der diagrams [|122[|123 ll. We can see some improvement of our results compared to naive mean-field 
approximations and other analytic approaches. 

Our e expansion predicts the behavior of the thermodynamic quantities near the unitarity limit. 
From Eqs. (13.441) . (13.451) . and (|3.42l) . the extrapolations to the three spatial dimensions 6=1 lead to 

' ^4-, (3.78) 



sfN d + 2 8 ep 



" 'f!l(5 l„£i|, (3.79) 



SpN d + 2 4eF\2 Sp 
ii«^_if^(3_lnf^l (3.80) 



£p 4 £p \ £p ; 

The pressure, energy density, and chemical potential at fixed density are decreasing function of the 
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binding energy s\, with the slope shown above. In Chap.|7l we will make a discussion to improve the 
extrapolation of the series expansion by imposing the exact result at J = 2 as a boundary condition. 
We can also try to determine the location of the splitting point. At this point, 

- = ^^^+0{e)]^2 atJ = 3 (3.81) 

and comparing with Eq. (12.141) . one finds that at the splitting point akp ~ 1. Since this result is 
known only to leading order in e, one must be cautious with the numerical value. However, certain 
qualitative features are probably correct: the splitting point is located on the BEC side of the unitarity 
limit (a > 0) and that the chemical potential fi/sp ~ -0.5 is negative at this point. 



Chapter 4 

Phase structure of polarized Fermi gas 



4.1 Background and proposed phase diagram 

Here we apply the e expansion developed in the previous Chapter to the unitary Fermi gas with 
unequal densities of two components (polarization). We put special emphasis on investigating the 
phase structure of the polarized Fermi gas in the unitary regime, which has a direct relation with the 
recent measurements in atomic traps [|^H3lM1,l45ll46l . Furthermore, it may be possible in future to 
realize the Feshbach resonances between two different species of fermionic atoms to study a Fermi gas 
with finite mass difference between different fermion species. Such asymmetric systems of fermions 
with density and mass imbalances will be also interesting as a prototype of high density quark matter 
in the core of neutron stars, where the density and mass asymmetries exist among different quark 
flavors [|62l l63l l64l l65l l66l . Extension of the e expansion to the polarized Fermi gas with unequal 
masses will be studied in the next Chapter and we concentrate on the equal mass case here. 

Reliable facts on the phase structure of the polarized Fermi gas is known only in the weak-coupling 
BCS and BEC limits. If the polarization chemical potential H = {p.^ - fUi)/2 is sufficiently small, the 
ground state is the unpolarized superfluid state both in the BCS and BEC limits and these two limits 
are considered to be smoothly connected as a function of 1/akp (phase I in Fig. 14.11) . When H is 
increased, we will have different situations in the two limits. In the BCS limit, the BCS superfluid 
state becomes unstable toward to the Fulde-Ferrell-Larkin-Ovchinikov (FFLO) states where Cooper 
pairs form with a nonzero momentum and the superfluid order parameter varies in space [|124[ I125II 
(phase IV). For the sufficiently large H, the Cooper paring is unfavorable and the system goes to a 
polarized normal Fermi gas (phase II). On the other hand, when H is increased in the BEC limit, 
unbound fermions are created on the top of the molecular BEC ground state. The system becomes 
a homogeneous mixture of the polarized fermions with the condensed molecules (phase III). For the 
sufficiently large H, all bound molecules disappear and the ground state is a fully polarized normal 
Fermi gas (phase II). Since the FFLO phase where the rotational symmetry is spontaneously broken 
does not appear in the BEC limit, such a phase should terminates somewhere between the BCS and 
BEC limits. 

From these knowledge in the two limits and assuming there is a phase transition line between 
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Figure 4.1: Proposed phase diagram in the plane of the diluteness parameter k = -l/ak^ and the po- 
larization chemical potential rj = H/Ah^q by Son and Stephanov II112L I is the unpolarized BEC/BCS 
phase, II is the normal phase, III is the gapless superfluid phase, and IV is a region of Fulde-Ferrell- 
Larkin-Ovchinikov phases ^1241 11251 . The dashed line divides phases II and III into regions with 
one (on the left) and two (on the right) Fermi surfaces. Region IV must be divided into phases with 
different patterns of breaking of the rotational symmetry (not shown). 



the phases III and IV instead of I and II, Son and Stephanov proposed a global phase diagram of the 
polarized Fermi gas in the BCS-BEC crossover as shown in Fig. 14.11 [ill2il . The point S where the 
FFLO phase terminates is called the splitting point and the phase structure around it is studied based 
on the effective field theory |112[| . The purpose of this Chapter is to study the phase structure of the 
polarized Fermi gas in the unitary regime including the splitting point from the microscopic point of 
view using the e expansion. 



4.2 Effective potential at finite superfluid velocity 

In order to take into account the possibility of the FFLO state in the phase diagram, we generalize 
our formalism to allow a spatially varying condensate where {(f){x)) = e^'"™'^0o with v^; being the 
superfluid velocity. The factor g^fnws a; jj^ fj-ont of 0o in the Lagrangian density can be absorbed by 
making the Galilean transformation on the fermion field as i/^o-(x) e""^'''"i/'o-(x) and the boson field 
as (fix) e^""'"'''^(p(x). Accordingly, the fermion propagator in Eq. (|3.7|) is modified as 



G(po,p) 



1 



(po + H -p- v^f - (Sp + Smv,)^ - (Pl + iS 



X 



Po + H + 

-00 



-^0 



(4.1) 
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while the boson propagator in Eq. (|3.8I) as 

D{po,p) = lpo-^-^ + iS] . (4.2) 



2 

The value of the superfluid velocity is determined from minimizing the dependent part of the 
effective potential Vh(Vs). As far as the polarization H is sufficiently small compared to the energy 
gap A of the fermion quasiparticle, VniVf^) to the leading order in e is given by the fermion one-loop 
diagram with the propagator in Eq. (|4.1l) : 

Vh<a{v,) = - J yjisp + s,^^f + 02. (4.3) 

Since it will turn out that Smv, is 0(6^'), we can expand the integrand in terms of £mvJ(po to lead to 

Vh<.(v.) = - J - — i^j • (4-4) 

Using the definition of the fermion number density A'^ in Eq. (13.401) . this part can be rewritten as 
Vh<a(v^) = Nsniv,, which represents the energy cost due to the presence of the superfluid flow. 

If H - p ■ reaches the bottom of the fermion quasiparticle spectrum A + {Sp + £,„v, - s^f jlcpQ, 
Vuiv^) receives an additional contribution from the filled fermion quasiparticles. Since s,,,^^ ~ e^^ 
is small, we can neglect it in the quasiparticle energy. Then the dependent part of the effective 
potential is given by 

/dp 
-— (H-p v,- coFip)X , (4.5) 
(Inf 

where we have introduced a notation (x)> = x^' 9{x) and co^ip) is the fermion quasiparticle spectrum 
^f(p) = a + (Sp - eo)2/20o derived in Eq. (|3.53l) . The effective potential VniVf^) to the leading order 
in 6 has the same form as that studied based on the effective field theory in 01121 . 

4.3 Critical polarizations 

Now we evaluate the effective potential Vh(Vs) at J = 4 as a function of Vs and H. Changing the 
integration variables to z = £p/0o and w = p ■ v^, we have 

2n 

X 



/m0() 

VHiv,) = i^ (4.6) 



^ <iz J" z Vl - |// - A - ^^—^^00 - 2 yJz(poSmv,w 



where zo = eo/0o ~ We can approximate -^/z0oe^ by ^jzo(poSmv, because the difference will be 
0{6^) and negligible compared to itself ^jzo^o^mv, ~ f^- Then the integration over z leads to 



VM = (— ) 00 — - -zo (— zo) Vl 



2 ^JSmv,(p0ZQ 



(4.7) 
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By introducing the dimensionless variables x and h as 



= xH^e]\', and ^LA = 2h(^e]\t„ (4.8) 



00 / 00 \37r 

we can rewrite the effective potential V//(t)s) = ^^(jc) in the simple form as 



(4.9) 



Now we see that if the superfluid velocity exists, Smv^ ~ e^^ and the dependent part of the effective 
potential is 0(6^°). 

Numerical studies on the effective potential Vh(x) as a function of x show that there exists a region 
of h, hi < h < h2, where Vhix) has its minimum at finite x. These two critical values are numerically 
given by 

hi = -0.00366 and h2 = 0.0275. (4.10) 



Correspondingly, we obtain the critical polarizations normalized by the energy gap as 
and 



^ = 1 - 0.0843 6^ 
A 



^ = l+0.63462f?f. (4.12) 



A \A/ 

Here we have replaced 0o by A because they only differ by 0(6) [Eq. (13.541) 1 . The region for the phase 
with spatially varying condensate is H2 - Hi ~ e^, where the superfluid velocity -Ug is finite at the 
ground state. 

As the polarization increases further H > H2, the superfluid velocity disappears. If H < 0)^(0), 
fermion quasiparticles which have momentum co^ip) < H are filled and hence there exist two Fermi 
surfaces, while there is only one Fermi surface for H > ojpiO) [see Fig. 13. 5L Therefore, from the 
quasiparticle spectrum derived in Eq. (13.531) . the polarization for the disappearance of the inner Fermi 
surface Hi, is given by 

^ = ^*(2) = i,i(^t (4.13) 

A A 2IA/ ^ ^ 

Here the location of minimum in the fermion quasiparticle spectrum so is related with the binding 
energy Sb near the unitarity limit via eo = eA- eb/2 [Eq. (|3. 581) 1. The critical polarizations Hi /A and 
H2/A as functions of Si,/eA are illustrated in Fig. 14.21 



4.4 Phase transition to normal Fermi gas 

Next we turn to the phase transition to the normal Fermi gas which occurs at // - A ~ 6. Since the 
region for the phase with spatially varying condensate is Hi 2 - A ~ e^, we can neglect the superfluid 
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velocity here. We can also neglect the possibility to have two Fermi surfaces where Ht, - Is. ~ e^. 
Then the contribution of the polarized quasiparticles to the effective potential is evaluated as 



(H-A)l (mcp,yi^ 



200 



In 



(4.14) 



where we have neglected the higher order corrections due to the shift of the location of minimum in 
the dispersion relation sq ~ e. According to the modification of the effective potential in Eq. (|3.36l) to 
^effC^^o) + Vh, the solution of the gap equation in Eq. (|3.38l) becomes 0o ^ 0o + 4>h, where 



(Ph = -(H - A)> + 



(.H-A)l 
200 



(4.15) 



Then the pressure of polarized Fermi gas in the superfluid state is given hy P = -Veff(0o + 0//) - Vh, 
which results in 



p=^ 
6 



y + ln2^ 



The polarization dependent part of the pressure Ph ~ is small and negligible in the region consid- 
ered here H - A ^ e. 

The pressure of the superfluid state should be compared to that of the normal state with the same 
chemical potentials. Since the phase transition to the normal state happens at // ~ 0o » yu ~ e(po, 
only one-component of fermions exists in the normal state. Therefore, the interaction is completely 
suppressed in the normal Fermi gas and its pressure Pn is simply given by that of a free Fermi gas 
with a single component: 



(H + ju)2+' 

r(f + 2) 



m \2 
2^1 



(4.17) 



The phase transition of the superfluid state to the normal state occurs at H = where the two 
pressures coincide P = Pn. From Eqs. (|4.16l) and (14.171) . the critical polarization satisfies the 
following equation, 



Sh ^ 2 + In 2 

1 + Ce — 6 

400 6 



00 



1 , /- 

4A 



2 + ln2 



-6 + (81n3- 121n2)6 



(4.18) 



where we have substituted the relation between the condensate 0o and the energy gap A at zero polar- 
ization in Eq. (13.541) . Defining a number cr ~ 0(1) by 



2 + ln2 

cr = C+ — (81n3 - 121n2) « 0.12197. 



(4.19) 
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the critical polarization normalized by the energy gap at zero polarization is written as 

f = 1 - ecr + ^ + O(e'). (4.20) 

If the binding energy is large enough Sh/eA > 4cr = 0.488, the superfluid state remains above H = A. 
The superfluid state at A < H < He is referred to as a gapless superfluid state. The fermion number 
densities of two diff"erence components are asymmetric there and the fermionic excitation does not 
have the energy gap. 

In particular, at the unitarity limit where e^l A = 0, the critical polarization is given by 



A 

At the splitting point where Sb/'^^ = 26, it is 



= 1 -6(T= 1 -0.122 e, (4.21) 

UL 



A 



= l-eo-+^ = l+ 0.378 6. (4.22) 

sp 2 



The extrapolations to three spatial dimensions 6=1 give the critical polarizations at the two typical 
points as //c/A|ul = 0.878 and //c/A|sp = 1.378. At unitarity, hence, there is no gapless superfluid 
phase. On the other hand, near the splitting point the normal phase is not competitive compared to the 
gapless phases. The phase boundary between the superfluid and normal phases HdAas a function of 
Eb/eA is illustrated in Fig. l4.2[ 



4.5 Phase structure near the unitarity Hmit 

The schematic phase diagram of the polarized Fermi gas in the unitary regime is shown in Fig. 14.21 
in the plane of H and s^je for the fixed energy gap A at zero polarization. The critical polarization 
Hc/A divides the phase diagram into two regions; the superfluid phase aiH < (I and III) and the 
normal phase aiH > (II). The Fermi gas in the normal phase is fully polarized near the unitarity 
limit because H :s> p. The phase transition of the superfluid state to the normal state is of the first 
order, because of the discontinuity in the fermion number density which is 0{e'^) in the superfluid 
phase while 0(1) in the normal phase. 

The superfluid phase can be divided further into two regions; the gapped superfluid phase atH < A 
(I) and the gapless superfluid phase at A < H < (III). At the BEC side of the splitting point (SP) 
where Sb/eA > 2, the phase transition from the gapped phase to the gapless phase is of the second 
order because a discontinuity appears in the second derivative of the pressure d^P/dH^ ~ e6{H - A) 
[see Eq. (14.161) 1. On the other hand, at the BCS side of the splitting point where 

4a- < ^ < 2, (4.23) 

6A 

there exists the superfluid phase with spatially varying condensate (IV) at < H < H2 between the 
gapped and gapless phases. This phase appears only in the narrow region where Hj- Hi ~ 6^. The 
phase transitions at H = Hi, H2 are of the first order. 
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Figure 4.2: Schematic phase diagram of the polarized Fermi gas near the unitarity limit from the e 
expansion. The phase diagram can be divided into four phases, I: the gapped superfluid phase, II: 
the polarized normal phase. III: the gapless superfluid phase, and IV: the phase with spatially varying 
condensate. The inset is the magnification of the region around the splitting point (SP). The phase IV 
appears in the narrow region represented by the thick line between the phases I and III. 

In actual experiments using fermionic atoms, unitary Fermi gases are trapped in the optical po- 
tential V(r) |l42l|43l|44l|45l|46]|. In such cases, the polarization H and the scattering length or e\, are 
constant, while the efl'ective chemical potential yu - V{r) decrease from the center to the peripheral 
of the trapping potential. For the purpose of comparison of our results with the experiments on the 
polarized Fermi gases, it is convenient to visualize the phase diagram in the plane of H and fi for the 
fixed Sb- 

In Fig. 14.31 eHc/sb and eA/sb are plotted as functions of /u/sb for the fixed binding energy > 0. 
The superfluid phase and the normal phase are separated by the line of the critical polarization H = 
while the gapped and gapless superfluid phases are separated by the line of the energy gap at zero 
polarization H = A. The intersection of two lines. He = A, is located at (yu/sb, e/Z/Sb) = ((1 - 
4cr)/8cr, l/4cr) = (0.524, 2.05), while the splitting point in the phase diagram is at (ji/sb, eH/sb) = 
(-1/4, 1/2). The horizontal dashed line from right to left tracks the chemical potential in the trapped 
Fermi gas from the center to the peripheral. If the polarization is small enough compared to the 
binding energy eH/sb < 2.05, there exists the gapless superfluid phase (III) between the gapped 
superfluid phase (I) and the normal phase (II). When the polarization is above the splitting point 
0.5 < eH/sb < 2.05, the phase with spatially varying condensate (thick line between I and III) will 
appear between the gapped and gapless superfluid phases. 

Assuming that the picture remains qualitatively valid in three dimensions, we thus found that the 
physics around the splitting point is the same as argued in Ref. |112| . However, due to the competition 
with the normal phase, the gapless phases disappear at some point, probably before the unitarity is 
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Figure 4.3: Schematic phase diagram of the polarized Fermi gas near the unitarity limit in the H-fi 
plane for the fixed binding energy > 0. The phases I, II, and III are same as in Fig. 14.21 The 
horizontal dashed line from right to left tracks the chemical potentials (pi, H) in trapped Fermi gases 
from the center to the peripheral. 

reached. The point where this happens can be estimated from our calculations to be 



The gapless phase with spatially varying condensate exists in a finite range 0.494 < {ak^y^ < 1- Thus 
the e expansion does not support the hypothesis of Ref. HI 1211 that the region of the gapless state with 
spatially varying condensate is connected with the FFLO region on the BCS side (see Fig. 14. II) . 



1 = V2^ = 0.494. 



(4.24) 



Chapter 5 

Fermions with unequal masses 



5.1 Two-body scattering in vacuum 

Now we generalize our discussion on the unitary Fermi gas in the e expansion to the fermions with 
unequal masses mi, especially to study the phase structure of unitary Fermi gas with unequal 
densities and masses. As we mentioned in the previous Chapter, such asymmetric systems of fermions 
with density and mass imbalances will be interesting as a prototype of high density quark matter in the 
core of neutron stars, where the density and mass asymmetries exist among different quark flavors ll62l 
l63l [64l [65l l66l . So far the phase structure of the two-component Fermi gas with density and mass 
imbalances in the BCS-BEC crossover has been studied within the mean-field approximation [|126[ 
I127[|128lll29ll . which is not a controllable approximation in the strong-coupling unitary regime. 

First we reconsider the two-fermion scattering in vacuum studied in Chap. [2l As a result of 
the summation of the geometric series of bubble diagrams, the inverse T-matrix of the two-body 
scattering is given by 

The limit of infinite scattering length corresponds to 1 /co ^ as before. Introducing the total mass 
M = nif + Ml and the reduced mass = m-|-m|/(m-|- -I- m^), the inverse T-matrix at d = 4 - e spatial 
dimensions to the leading order in e can be written as 

r(p.,p)-. = i-i(^|)^(p„-|l. ,•*) + .... (5.2) 

This result suggests that the boson propagator is given by 

D(po,p) = Lo-^ + is\ (5.3) 



with the fermion-boson efi'ective coupling 



(5.4) 
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The binding energy £b > is obtained as the location of the pole of the T-matrix at zero external 
momentum: T(-£b, OT^ = 0. From Eq. (15.11) . we see cq is related with eb via the following equation 



(5.5) 



Near the four spatial dimensions, the relationships between eb and cq to the next-to-leading order in 
6 = 4 - J is given by 



1 

Co 



'2e\7T 



i + 2(i-r) 



~2n 



(5.6) 



5.2 6 expansion with unequal masses 

Now we consider the finite density system. The formalism of the e expansion developed in the equal 
mass case in Chap. Oholds just by replacing the boson propagator D{p) by Eq. (15.31) and the fermion- 
boson effective coupling g by 

{In^efl^ (m,<Po\^l' ^^^^ 



(5.8) 
(5.9) 
(5.10) 



Then the Lagrangian density under consideration is given by 



= g^V+¥^ + g¥V_¥/ + /z^FVa^ + \2iu-—\<p*<p- —if - —if* 



Co 



Co 



Co 



where fXa- = IJ^±H and l/lmu- = (1 ± K)/4my. We introduced the dimensionless parameter -1 < k < 1 
defined by 

K 



m-f + mi 

which measures the mass difference between two fermions. 

The Lagrangian density Xo generates the fermion propagator given by 



(5.11) 



G(pq,p) 



1 



(Po+H - KSpf -E1 + id 



PO + H - KSp + £p 
-00 



-(po 

Po+H -KSp- £. 



(5.12) 



where Sp = p^/Ain,- and Ep = + cp^, and the boson propagator given by 



(5.13) 



We can use the same Feynman rule and the same power counting rule of e as the equal mass case 
developed in Sec. 13. 2[ as far as the mass difference is not so large. When the mass difference is 
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as large as \k\ - I - 0(1/6), the power counting rule of € breaks down because of the existence 
of another large parameter mheavy/'^iight ~ 1/f- Accordingly, the fermion-boson scattering leads to 
the three-body bound states (Efimov states) in vacuum when the mass ratio is given by mheavy/"^iight = 
4/6+0(1) [|134[|135lll36L In this Chapter, we will only consider the case where the mass difference is 
at most /c ~ 0(1) and free from the Efimov states. First we will determine the pressure and the fermion 
quasiparticle spectrum to the leading and next-to-leading orders in e at the unpolarized ground state 
where H = 0. Then we investigate the phase structure of polarized Fermi gas with the small mass 
difference k ~ e. 



5.3 Effective potential and pressure 

The computation of the effective potential with finite k is straightforward according to the systematic 
expansion over e. Since KSp < Ep for any momentum p, k does not appear at the one-loop level. 
Therefore, the effective potential from the tree and one-loop diagrams depends only on the reduced 
mass rtij- as 

Vo{(f>o) + Vi(0o) 



,2 



Co 



00 L ^ 7-3(y + ln2) 1 l-2(y-ln2) ^ 



d/2 



3 [ 6 j 6 

The effective potential at two-loop level is given by the same diagram in Fig. l3.3l as 

V2(0o) = / J Gnip)G22iq)Dip - q) 

-2 Cdpdq {Ep-ep){Eq-Sq) (5.15) 



r dpdq 
4 J (Inf'' 



EpEq ^Ep -|- Eq -\- 2 ^p~q ^^p ^^(^ 



Changing the integration variables io x = £p/4>o, y = and cos 6 = p ■ q, the integral can be 

expressed by 

,2 



m 

V2((f>o) = -e 



rhv^o r, r. r,, . 2. [m-x][f(y)-y] 

— — dx \ dy \ dOxy sin 9 f r 

n I n Jo Jo Jo f (x)f{y)[g,{x,y)- (I -K^)^ cose] 

(5.16) 

with f(x) = Vx^TT and gi,(x, y) = f(x) + f(y) + ^-^(x + y) - k(x - y). The integration over 6 can be 
performed analytically to lead to 

/mi.0o'^ 



V2(0o) = -C,6^^ I <^o (5.17) 



where is a function of k given by a two-dimensional integral 



^0 
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is symmetric under k ^ -k and numerically we find C^^ is an increasing function of a: at < a: < 1 . 
In the equal mass limit /<• = 0, we reproduce Co ~ 0.14424 in Eq. (13.351) . 

Consequently, we obtain the effective potential Veff = Vq + + ^2 up to the next-to-leading order 
in 6, 



Veff(<^0) = — + 
Co 



7-3(7 + ln2) 



1 -2(7-ln2) 



n 



d/2 



(5.19) 



Comparing Veff(0o) with the effective potential in the equal mass case in Eq. (13.361) . we see the solution 
of the gap equation (po is simply given by replacing the fermion mass m by 2m^ and the constant C by 
C^. Using the relation of cq with the binding energy of boson Sb in Eq. (|5.5I) . we obtain the condensate 
as 



= ^[1 +(3C,- 1 +ln2)6] + — 
e e 



l+|3C,-i+ln2-^ln^|6 



Then the pressure up to the next-to-leading order in e is given by 



P = -Kff(0o) = ^ 



,17 y + ln2 



3£b 

400 



n 



d/2 



(5.20) 



(5.21) 



When the mass difference k is as small as e which is the case we will concentrate on later, we can 
neglect the k dependence in C^, to give C ~ 0.14424. In this case, the pressure of the superfluid state 
up to the order 0{e) depends only on and not on k at all. 



5.4 Self-energy and dispersion relation 

Let us turn to the computation of the fermion quasiparticle spectrum with finite k. To the leading 
order in e, the dispersion relation of fermion quasiparticles are given by 

<^f(p) = Ep± KSp. (5.22) 

The lighter fermion has the energy gap A = 0o at Sp = 0, while the heavier fermion has the energy 
gap A - (pQ at Sp = \k\ 4>oI Vl - K^. The next-to-leading order corrections to the dispersion 

relation come from the jj. insertion to the fermion propagator and the one-loop self-energy diagrams, 
-il.{p), depicted in Fig. 13. 4[ The one-loop diagrams in Fig. 13.41 give corrections to the diagonal 
elements of the fermion self-energy, which is evaluated as 



£ii(p) = ig' j G22(k)Dip - k) 



2 J {2ny 



yl+l 

Ek - Sk 



(5.23) 



Ek [Ek + KSk + ^£k~p - Po] 
and 

^22(p) = eJ^Gn(k)D(k-p) 

dk Ek - Sk ^^•^'^^ 



2 J (2n 



y Ek [Ek - KSk + ^Sk-p + Po] 
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From now on, we concentrate on the case where the mass difference is as small as k ~ e, where 
we will see interesting physics in the phase diagram. In this case, we can neglect the k dependence 
in Ell and £22 to the order 0(6) because Z is already small by the factor ~ e. Therefore, the mod- 
ification of the fermion quasiparticle spectrum due to the finite k comes from Eq. (15.221) . Following 
the calculations in Sec. 13.51 we find that the dispersion relation of the fermion quasiparticle around its 
minimum has the following form 

coAp) - a + ^ 7" — - J(Sp ± #0 - eoY + A2. (5.25) 

200 ^ 

Here A is the energy gap of the fermion quasiparticle defined in Eq. (|3.54l) and sq is given by 
Eq. (13.551) . Note that the energy gap is not affected by k up to the order 0(6). The minimum of 
the dispersion curve is located at a nonzero value of momentum |p| satisfying Sp = sq - |/c|0o for the 
lighter fermion and Sp = sq + |/c|0o for the heavier fermion when those are positive. 



5.5 Critical polarizations and phase diagram 

Since the result should be symmetric under H -H and k -k, we can choose H > without 
losing generality. Accordingly, a- > (< 0) corresponds to the system where the majority is the 
lighter (heavier) fermions and an increase of k means a decrease of major fermion's mass. When H is 
increased, the phase transition to the normal Fermi gas occurs at // - A ~ 6. The critical polarization 
He is given where the pressure of the superfluid state coincides with that of the normal state with 
the same chemical potentials. The pressure of the superfluid state is given by Eq. (15.211) because the 
contribution of the polarized quasiparticles to the pressure is 0{e^) and negligible [see Eq. (|4. 161) 1. 
On the other hand, since // ~ 00 » ~ 60o, the pressure of the normal state is given by that of the 
fully polarized free Fermi gas 

The condition P = for the phase transition to the normal state gives the critical polarization H = 
normalized by the energy gap at zero polarization A in Eq. (|3.54l) as 

He Sh 2 T 

-f = l-ecT+^ + -K+0{e\ (5.27) 

The number a ~ 0.12197 is defined in Eq. (14.191) . When > A, the gapless superfluid state appears 
at A < // < He- In particular, in the unitarity limit £5 = the gapless superfluid state is possible when 
the majority is the lighter fermions and the mass difference is as large as K/e > 3cr/2 = 0.183. The 
phase boundary between the superfluid and normal phases He/A. in the unitarity limit as a function of 
K/e is illustrated in Fig. 15.11 

The minimum of the dispersion curve for the major fermions is located at a nonzero value of 
momentum \p\ satisfying 

so-KcPo = {e-K)A-^ + 0{e^). (5.28) 
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H/A 




K 

1 1 1 1 >- — 

-0.5 0.5 1 ^ 

Figure 5.1: Schematic phase diagram of the polarized Fermi gas with unequal masses in the unitarity 
limit from the e expansion, /c > (< 0) corresponds to the system where the majority is the lighter 
(heavier) fermions and an increase of k means a decrease of major fermion's mass. The phase diagram 
can be divided into four phases, I: the gapped superfluid phase, II: the polarized normal phase. III: 
the gapless superfluid phase, and IV: the phase with spatially varying condensate. The inset is the 
magnification of the region around the splitting point (SP). The phase IV appears in the narrow region 
represented by the thick line between the phases I and III. 



Therefore, we find the splitting point where the minimum of the dispersion curve sits exactly at zero 
momentum exists at 

— + K = e. (5.29) 
2A ^ 

The phase structure around the splitting point is determined by the same discussion given in Chap. E] 
just by replacing so by eo - k^Po- There exists the phase with spatially varying condensate (or finite 
superfluid velocity) at Hi < H < H2, where 



^ = 1 - 0.0843 6^ 
A 



j-x)', (5.30) 



and 



^ = 1 +0.634 6^ (I -/c)'. (5.31) 
Furthermore, the polarization for the disappearance of the inner Fermi surface Hj, is given by 

«3 ".(0)^, 



A A 2\A 

If A < H < Ht,, fermion quasiparticles which have momentum co^ip) < H are filled and hence there 
exist two Fermi surfaces, while there is only one Fermi surface for H^ < H. Since so-K(f>o = (e - k) A 
in the unitarity limit Sb = 0, the splitting point is located at k = €. The critical polarization for the 
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Figure 5.2: Phase diagram of the polarized Fermi gas in the plane of binding energy and mass 
difference k from the e expansion. Two solid lines show = A in Eq. (15.341) and the location of 
the splitting point (SP) in Eq. (15.291) . The intersection of these two lines is given by Eq. (15.351) . The 
phases II, III, and IV represents the ground state slightly above H = A. The phase diagrams in the 
equal mass limit k = and in the unitarity limit Sb = with //-axis are already illustrated in Fig. 14.21 



and Fig. l5.1[ respectively. 

phase transition to the normal state at this point is Hc/A = 1 + 0.5456 > 0. Therefore, the phase with 
spatially varying condensate and the phase with two Fermi surfaces are stably exist in the unitarity 
limit with finite mass difference given by 



The critical polarizations Hi /A and //2/A in the unitarity limit as functions of k/€ are illustrated in 



Finally, we discuss the whole phase structure of the polarized Fermi gas with unequal masses near 
the unitarity limit. Hdsb, k) in Eq. (15.271) is a function of the binding energy Sb and the mass difference 
K. Hence H = forms a single plane in the three-dimensional phase diagram in the space of Sb/A, k, 
and HI A, which separates the superfluid phase at H < He from the polarized normal phase dXH > H^. 
The splitting point given by Eq. (15.291) with H = A forms a "line" in such a three-dimensional phase 
diagram. The superfluid phase with spatially varying condensate emerges from the SP line, which 
occupies a finite region in the phase diagram between the gapped superfluid phase at // < //i and the 
gapless superfluid phase at H2 < H < H^. 

Fig. [521 shows a part of such a three-dimensional phase diagram in the space of s^jA, k, and HI A, 
sliced slightly above H = A. When He < A for a given set of Sb and k (lower half of the He = A line 




2 6 



(5.33) 



Fig.O 
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Figure 5.3: Schematic phase diagram of the polarized Fermi gas in the H-s\, plane near the disappear- 
ance of the splitting point at /c = -2.0 e and -2.5 e. The phases I, II, III, and VI are same as in Fig. 15. II 
If K > -2.27 6 (left panel), the splitting point and the phase with spatially varying condensate stably 
exist, while those become unstable when k < -2.27 e (right panel). 



in Fig. 15.21) . the superfluid state is unstable at // = A and there exists the phase transition from the 
gapped superfluid state to the polarized normal state (II). On the other hand, when > A (upper half 
of the He = A. line in Fig. 15.21) . the superfluid is stable even at // > A and the gapless superfluid state 
(III) appears there. The boundary He = A dividing the phase diagram in the plane of eb and k into 
these two regions is given from Eq. (|5.27l) by 

^ + = (5.34) 

As the mass of major fermions is increased (or k is decreased), the boundary shifts to the larger value 
of £b because the pressure of the normal state increases. Consequently, the polarized normal state 
occupies the wider domain near the unitarity limit. 

The region where H^. > A can be further divided into two regions according to whether the min- 
imum of the dispersion curve is located at a nonzero value of momentum sq - K(pQ > or not. If 
eo - K^(pQ > (below the SP line in Fig. 15.21) . there exist the superfluid state with spatially varying 
condensate (IV) between the normal state and the gapless superfluid state. The boundary sq = K(f>o 
represents the line of the splitting point given by Eq. (15. 291) . As the mass of major fermions is in- 
creased, the region for the superfluid phase with spatially varying condensate shrinks. Eventually, 
such a phase disappears at the following point where the two hnes He = A and £o = #0 merges 

^ =4(2- 3o-)e = 6.546, 

^ (5.35) 
K = -3 {I - 2(t) e = -2.27 e. 
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Therefore, if the mass difference is as large as a: < -2.27 e with major heavier fermions, the splitting 
point becomes unstable towards the polarized normal state and the superfluid phase with spatially 
varying condensate never appears for any value of the binding energy s^,. The schematic phase dia- 
gram in the plane of Sb and H near the disappearance of the splitting point is illustrated in Fig. I5.3l by 
taking the typical values at k = -2.0 e and -2.5 e. 

5.6 Summary on the phase structure of polarized Fermi gas 

In Chapters |5] and |5l the phase structure of the polarized Fermi gas with equal and unequal fermion 
masses has been studied in the unitary regime based on the e expansion around four spatial dimen- 
sions. Although our results are valid only near four spatial dimensions where e <^ 1, we can draw 
the following conclusions assuming that the picture remains qualitatively valid in three dimensions. 
At unitarity in the equal mass limit, there is a first-order phase transition from the unpolarized su- 
perfluid state to the fully polarized normal state. On the BEC side of the unitarity point, the gapless 
superfluid phase and the superfluid phase with spatially varying condensate stably exist between the 
gapped superfluid phase and the polarized normal phase in a certain range of the binding energy and 
the mass difference. In the equal mass limit (Fig. 14.21) . our study gives a microscopic foundation to the 
phase structure around the splitting point, which has been proposed on the BEC side of the unitarity 
point using the effective field theory [1121. Moreover, we found the gapless phases terminate at some 
point before the unitarity is reached due to the competition with the polarized normal state. The naive 
extrapolation of our result (14.231) with the use of Eq. (12.141) estimates the range of the superfluid phase 
with spatially varying condensate to be 0.494 < (afcp) ' ^ 1 at three dimensions 6^1. 

We also found the splitting point and the same phase structure around it in the unitarity limit but 
with finite mass difference between two fermion species when the majority is the lighter fermions 
(/<■ > in Fig. 15.11) . The range of the superfluid phase with spatially varying condensate can be 
estimated from Eq. (|5.33l) to be 1.37 < mminonty/'^majority < 3 at 6 ^ 1. Such splitting points form a 
smooth line on the H = A plane of the three-dimensional phase diagram in the space of Sb/A, k, and 
H/A. It gets away from the unitarity with increasing the mass of major fermions (or decreasing k in 
Fig. 15.21) . Eventually the splitting point becomes unstable due to the competition with the polarized 
normal state at the point in Eq. (15.351) . which we estimate to be (ak^y^ « 1.81 and mmajority/''^mmority ~ 
5.54. Accordingly, the region for the superfluid phase with spatially varying condensate shrinks and 
eventually disappears. The superfluid state with finite superfluid velocity is never realized for any 
value of the binding energy Sb for a sufficiently large mass difference where the majority is the heavier 
fermions mmajority/'^minority ^ 5.54. While our quantitative results are not reliable at three dimensions, 
it is interesting to point out that our estimates on the mass ratios are well below the critical value at 
d = 3, mheavy/'J^iight = 13.6, whcrc the instability occurs due to the Efimov effect [11341 11351 1136II . 
Further study will be worthwhile to confirm these possibilities. 
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Chapter 6 

Expansion around two spatial dimensions 



6.1 Lagrangian and Feynman rules 

In this Chapter, we formulate the systematic expansion for the unitary Fermi gas around two spatial 
dimensions in a similar way as we have done for the e = 4 - d expansion. Here we start with the 
Lagrangian density given in Eq. (|3.2I) limited to the unpolarized Fermi gas in the unitarity limit where 
// = Oand 1/co = 0: 

X = (idt + + JU0-3) »F + ^V+^F^ + *FV_¥(^*. (6.1) 



2m 

Then we expand the field around its vacuum expectation value 00 as 

= 00+^^, ^ = l2^) ' ^^-^^ 

where the effective coupling g ~ e^^^ in Eq. (|2.10l) was introduced. The extra factor (m/j./lny'^''^ was 
chosen so that the product of fields (p*(p has the same dimension as the Lagrangian density §. 
Then we rewrite the Lagrangian density as a sum of three parts, X = Xo + Xi + X2, where 

Xo = ^idt + + IU0-3 + + cr-0()j ^ , (6.3) 

Xi = -(p*(p + ^*fV+¥^ + ^¥V^»P^*, (6.4) 

X2 = <^>. (6.5) 
The part Xo represents the gapped fermion quasiparticle, whose propagator is given by 



G(po,p)= ^ 



^Po + Sp- p -00 



-00 Po-ep + 



(6.6) 



'The choice of the extra factor is arbitrary, if it has the correct dimension, and does not affect the final results because 
the difference can be absorbed into the redefinition of the fluctuation field if. The particular choice of g in Eq. ( 16.21 ) will 
simplify the form of loop integrals in the intermediate steps. 
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propagators 



iG 



vertices 



tga+ 
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Figure 6.1: Feynman rules for the expansion around two spatial dimensions from the Lagrangian 
density in Eq. (I6.3|) . The first line gives propagators, while the second line gives vertices. 



with Ep = JiSp - nY + 00 being the usual gapped quasiparticle spectrum in the BCS theory. 



The second part Xi represents the interaction between fermions induced by the auxiliary field ip. 
The first term in Xi gives the propagator of the auxiliary field ip 



and the last two terms give vertices coupling two fermions and ip. If we did not have the part X2, we 
could integrate out the auxiliary fields ip and cp* to lead to 



which gives the contact interaction of fermions with the small coupling ~ 6 as depicted in Fig. 12.11 
The vertex in the third part X2 plays a role of a counter term so as to avoid double counting of a 
certain type of diagram which is already taken into Xi as we will see below. The Feynman rules 
corresponding to these Lagrangian densities are summarized in Fig. 16. 1[ 

6.2 Power counting rule of e 

We can construct the similar power counting rule of e as in the case of the expansion around four 
spatial dimensions. First, we make a prior assumption 0o//^ ~ e'^^'^, which will be checked later. 
Since e"'^^ is exponentially small compared to any powers of 6, we can neglect the contributions of 0o 
when we expand physical observables in powers of e. Then, since each pair of fermion and (p vertices 
brings a factor of e, the naive power of e for a given diagram is Ngjl, where Ng is the number of 
couplings g from Xi . 

However, this naive counting does not take into account the fact that there might be inverse powers 
of e that come from integrals which have logarithmic divergences dXd = 2. Each loop integral in the 
ultraviolet region behaves as 




D(po,p) = -1, 



(6.7) 



(6.8) 




(6.9) 
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+ = 0(e) 

a) (b) 




(c) 

Figure 6.2: Two apparent exceptions of naive power counting rule of 6, (a, c). The boson self-energy 
diagram (a) is combined with the vertex from X2 (b) to restore the naive e counting. The condition of 
disappearance of the tadpole diagram (c) gives the gap equation to determine the value of condensate 

00- 



while each fermion propagator behaves, at worst, as G{p) ~ p ^. Therefore, a given diagram may 
diverges as ~ with D being the superficial degree of divergence given by 

D = 4L-2Pf- (6.10) 

Here L is the number of loop integrals and Pp is the number of fermion propagators. Using the similar 
relations to Eq. (13.131) . 

L = Pp + PB-N-g + l, 

E, ^^-''^ 
A^, = Pf + Y = 2Pb+£b, 

the superficial degree of divergence is written in terms of the number of external fermion (auxiliary 
field) lines, ^fcb). as 

£) = 4-£f-2£b- (6.12) 

Therefore, the inverse power of e is possible only in diagrams which satisfies Ef + 2£'b ^ 4. Moreover 
from the analytic properties of G(p) in the ultraviolet region discussed in Eq. (13.161) . one can show 
that there are only two skeleton diagrams which have the 1/e singularity near two dimensions. They 
are one-loop diagrams of the boson self-energy rFigs. l6.2l' a)1. and the if tadpole diagram rFig. l6.2r c)1. 
The boson self-energy diagram in Fig. l6.2r a) is evaluated as 



(6.13) 



-■n.0') = -i^/(5^«n(*.f)G4-f) 

_ ._2 r dk 0(Sk+R - P) - - Sk-R) 



i.Po - + 2// -I- iS) 
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where we have neglected the contribution of 00. The integral over k is logarithmic divergent atd = 2 
and has a pole at 6 = 0. Thus it is 0(1) instead of 0{e) according to the naive counting. The residue 
at the pole can be computed as 

_ _^ r dk 1 
n.ip) J ^2ny 2sk - {po - \ep + 2^ + i5) ^ " " 

= 1 + 0(6), (6.14) 

which is cancelled out exactly by adding the vertex Hq = -1 in X2- Therefore the diagram of the type 
in Fig. I6.2r a) should be combined with the vertex from X2 in Fig. I6.2r b) to restore the naive e power 
counting result, 0{e). 

Similarly, the tadpole diagram in Fig. l6.2r c) also contains the 1 / e singularity. The requirement that 
this tadpole diagram should vanish by itself gives the gap equation to the leading order to determine 
the condensate (po as wee will see in the succeeding section. 



6.3 Gap equation 

The condensate (pQ as a function of the chemical potential jx is determined by the gap equation which 
is obtained by the condition of the disappearance of all tadpole diagrams. The leading contribution to 
the gap equation is the one-loop diagram drawn in Fig. l6.2r c). which is given by 



By changing the integration variable to z = Spljx, we obtain 



r(f) Jo ^ ^(z-ir + ih/f^y 



^, = ^^"^ Hdz ^''^ (6.16) 

' ^ m Jo ■ ■ 



The integration over z in the dimensional regularization can be performed to lead to 

(6.17) 



"1 = - ^ 



In — - - + 0{e') 



The first term In 2p/(pQ originates from the singularity around the Fermi surface as is well known as the 
Cooper instability, while the second term 1/6 is from the logarithmic singularity of the k integration 
at = 2 in Si. Solving the gap equation Si = 0, we obtain the condensate as (pa = 2fie~^^'^. Note that 
this result is equivalent to that obtained by the mean field BCS theory. 

It is known that the pre-exponential factor in the mean field result (po = 2fj. e'^^'^ is modified due 
to the effects of medium [11301 [TTlil . In the language of the tadpole diagrams, the corresponding 
modification to the gap equation comes from the three-loop diagram S3 depicted in Fig. 16.31 The 
diagram, which seems proportional to ~ 6^, gives the 0(1) correction to the gap equation. 



6.3. Gap equation 
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Figure 6.3: The tadpole diagram which gives the medium-efFect correction to the gap equation. 



Using the Feynman rules, the tadpole diagram in Fig. l6.3l is given by 

= -f J ^11^^^ ^22(P) Gn{p - k) G2M) G22(q + k). 



(6.18) 



Since the second term in the product, 



Gu(p)G22{p) = 



1 



+ 



(6.19) 



gives only the 0{e^) correction to the gap equation, we can neglect it for the current purpose. Then 
we perform the integration over po and to result in 



Cdkd 



dp dq <pQ 



Ep k() + Sp ic jj. ^ Ep + Sp ji 



(ko + Epy-E- 



X 



(Ep k^Y E^ j^ v-^u ' ^p/ ^p+k 

kp + Eq- £q+k + M Eg- Eg + p. 

{h + Egf - El^ ^ (ko - Egf - El^ 



(6.20) 



Because of the factor 1 /EpEg in the integrand, the integrations over p and q are dominated around the 
Fermi surface where ep(q) ~ /i and hence Ep^g^ ~ (pQ. Keeping only the dominant part in the integrand, 
we can write the integral as 

; dp dq 0o 1 



/dk dp dq 
(2;r)3^+i 



AEpEg (ko + - P){h) + Sg^k - P) 

Now the integration over ko can be performed easily to lead to 

dk 0{Sq+k - p)6{p- Sp-k) 



(6.21) 



^ ,g Cdpdq 1 r dk 
■^3 = -^2-0oJ ^^J^J J2n)" 



^q+k ^p-k 



(6.22) 



If we evaluated the k integration at (i = 3, we would obtain the static Lindhard function representing 
the medium- induced interaction f2T\. 

Here we shall evaluate E^atd = 2. Changing the integration variables to Sp, Sq, Sk, cos^p = k p, 
cosxg = k ■ q, and performing the integrations over Sp and Sq, we obtain 

2 

r dei. dy„ dv., 6(£„^u — u) 6(u — s^^u) 

(6.23) 



^q+k ^p-k 
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The range of the integrations over;^fp sindxq are from to n. Since cpof/J. oc e"'^*^, we have (2 lnyu/0o)^ ~ 
(2lef which cancels coming from the four vertex couplings g^. Finally the integrations can be 
performed as follows: 

dsk dXp dXq 6(£q+k -f^)9(M- £p-k) 



rdskdx 
J tt' 



Jo '"Jo ^ Jo ^ 



+ 

1 

2' 



2^JJIs^(cosXp + cosXq) 

J4ucos2y„ Jo ^ Jff/2 ^ 



(6.24) 



which gives S3 as 



m 



S3 ^ -ig(l)Q—. 

2n 



(6.25) 



Consequently, the gap equation Si + S3 = 0, which receives the 0(1) correction due to S3, is 
modified as 



ln^-i-l + C>(e) = 0. 
00 e 



(6.26) 



The solution of the gap equation becomes 



00 = 2/i exp 



-3-1+0(6) 

e 



^[1+0(6)]."^^ 

e 



(6.27) 



where the value of condensate is reduced by the factor e ~ 2.71828. The reduction of the pre- 
exponential factor due to the medium effects is known as the Gor'kov correction at J = 3 theo- 
ries ifBOlfBTH . 



6.4 Thermodynamic quantities 

The value of the effective potential Veff at its minimum determines the pressure P = -V^sicpo) at a 
given chemical potential //. Since the energy gain due to the superfluidity 0q ~ e~^^'^ is exponentially 
small compared to any power series of e, we can simply neglect the contributions of 0o to the pressure. 
To the next-to-leading order, the effective potential receives contribution from two vacuum diagrams 
drawn in Fig. 16. 4t fermion loops without and with an exchange of the auxiliary field. The one-loop 
diagram at = is 0(1) and given by 

Vm = -2 r -^OU - ep)> = -^J^ ^ -Pfree, (6.28) 

J W r(f + 2) 



6.5. Quasiparticle spectrum 
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0(1) 0(e) 
Figure 6.4: Vacuum diagrams contributing to the pressure up to next-to-leading order in e. 



which represents the contributions of free fermions to the pressure. The two-loop diagram at 0o = 
is 0(6), which represents the density-density correlation as 

dp dq 



_2 C dp dq 
J (27r)2^+2 



Gu{p)G22{q) 



/dp 



= -6- 



(6.29) 



r(f + if 



Thus we obtain the pressure up to the next-to-leading order in e as 

/' = (! +6)/'free. 



(6.30) 



Accordingly, the fermion number density is given by A'^ = dPjdp. = (1 -i- e)Nf^sc- The Fermi 
energy is obtained from the thermodynamics of free gas in d spatial dimensions as 



In 

Sf = — 
m 



= (1 +6)yU, 



which yields the universal parameter of the unitary Fermi gas from the e expansion as 

^=- = 1-^+0(6^). 



(6.31) 



(6.32) 



6.5 Quasiparticle spectrum 

To the leading order in e, the dispersion relation of the fermion quasiparticle is given by covip) = 
Ep = -sjisp - pY + (f>Q, which has the same form as that in the mean field BCS theory. There exist the 
next-to-leading order corrections to the fermion quasiparticle spectrum from the one-loop self-energy 
diagrams, -z£(p), depicted in Fig. 16.51 These corrections are only to the diagonal elements of the 
self-energy and each element is evaluated as 

dk 



^nip) 



-2 r 

_2 r dk 



G22(k) 



(6.33) 



^0(p- Sk) = -en 
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-ih = I \ + / \ 

> I ^ I ^ , I . L 



Figure 6.5: One-loop diagrams contributing to the fermion self-energy in order 0{e). 



and 



^ ^ ' (6.34) 



_y r dk _ 



To this order, the self-energy is momentum independent, which effectively shifts the chemical poten- 
tial due to the interaction with the other component of fermions. 

By solving the equation det[G '(a*, = in terms of to, the dispersion relation of the fermion 

quasiparticle up to the next-to-leading order is given by 

ojpip) = ^{Sp- n-enf + (pl- (6.35) 

The minimum of the dispersion curve is located at a nonzero value of momentum, \p\ = (2meo)^^^, 
where 

So = (1 + e))"- (6.36) 

The location of the minimum coincides with the Fermi energy in Eq. (|6.31l) . Eq = Sj^, in agreement 
with the Luttinger theorem [il32il . The energy gap A of the fermion quasiparticle is given by the 
condensate, 

A = 00 = — e'"'. (6.37) 
e 



6.6 Extrapolation to e=l 

Now we discuss the extrapolation of the expansion over e = d -2to the physical case at three spatial 
dimensions. In contradiction to the case of 6 = 4-J expansion, the coefficients of 0(e) corrections are 
not small. If we naively extrapolate the leading and next-to-leading order results for ^ in Eq. (|6.32|) . 
A in Eq. (|6.37l) . and sq in Eq. (16.361) to e = 1, we would have 

^«0, -« 0.271, — «2, (6.38) 

which are not as good as the extrapolations in the expansions over e in Eq. (13.771) . Thus, instead of 
naively extrapolating the e expansions to J = 3, we use them as boundary conditions to improve the 
series summation of the e expansions in Chap.|7l 



6. 7. NNLO correction for ^ 
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Figure 6.6: Vacuum diagram contributing to the pressure to the order e^. The counter vertex -iHo = i 
for each bubble diagram is understood implicitly. The other 0{e^) diagrams identically vanish. 



6.7 NNLO correction for f 



The order contribution to the effective potential or the pressure is solely given by the three-loop 
diagram depicted in Fig. l6.6[ which is written as 

dp ^(S+l -f^)-0(M- ep-|) 



(6.39) 



r)'' 2sp - (ko - + 2jj. + id) 

where po integrations in each bubble diagram are already performed. Note that +1 in the bracket 
comes from the counter vertex -Ho = 1 . We can see that the other 0(6^) diagrams identically vanish. 
The integration over ko in leads to 



^3(0) 



_2 r dkdp 



1 + 



^ J w 



dq -MSg-^ -f^) 



2s - 2s 



(6.40) 



Due to the ^-functions, the range of integral over s^ is limited to < Sfc < 4// and the range of integral 



over Sp is limited to < < Ap, where 



I cosxp 



Sk sm Xp 



with co^Xp = k ■ p. Similarly, the range of integral over Sg is limited to A,^ < Sq where 



(6.41) 



I cos^^ 



£k sm Xq 



with cos;^^ = k ■ q. Then the integration over Sg is performed to result in 



471- Jo Jo :t Jo n Jo 



y + In 



A - e 



(6.42) 



(6.43) 



Finally, introducing the dimensionless variable z = Sk/l^ and performing the integration over Sp, we 
obtain the following expression for V3, 



V,(0) = e 



17 



r + ^Jzj^ ^|A,lnA,-(A,-Ap)ln(A,-Ap)-Aj 



(6.44) 
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where Ap(^) = Ap(^)//i. The numerical integrations over z,;^'p, and;^^ give 

%(0) = e^^il+ 0.0568528] . (6.45) 
2n \2 I 

Therefore, combining this result with Eqs. (I6.28|) and (16.291) . we obtain the pressure up to the 
next-to-next-to-leading order (NNLO) in e as 

= Pfree [l + 6 - 0.3068528 6^] . (6.46) 

Then ^ up to the order is found to be given by 

^= — = [l + 6- 0.3068528 e^\^' = \-e+ 1.80685 + 0{e^). (6.47) 

Although the 0{e^) correction to the pressure is small, the NNLO correction to ^ turns out to be large 
because of the large 0{e) correction in the pressure. 



free 



1+6 + 



Chapter 7 



Matching of expansions around d = 4 and 

d = 2 



As we have mentioned previously, we shall match two expansions around d = 4 and d = 2, which are 
studied in Chapters [3] and[6]respectively, in order to extract results at J = 3. We use the results around 
two spatial dimensions as boundary conditions which should be satisfied by the series summations 
of the expansions over e = A - d. Because we do not yet have a precise knowledge on the large 
order behavior of the expansion around four spatial dimensions, we assume its Borel summability 
and employ Pade approximants. 

Let us demonstrate the matching of two expansions by taking ^ as an example. In Ref. Il96l . the 
linear interpolation between exact values atd = 2(^= 1) and J = 4 = 0) was discussed to yield 
^ = 0.5 at J = 3. Now we have series expansions around these two exact limits. The expansion of ^ 
in terms of 6 = 4 - J is obtained in Eq. (13.431) . Assuming the Borel summability of the e expansion, 
we write ^ as a function of e in the form of the Borel transformation, 

e^'^ / 6ln6 \ 

^(6) = — explg-^^jj^ dte-'B^(et), (7.1) 

where we factorized out the non-trivial dependence on e explicitly. B^{t) is the Borel transform of the 
power series in ^(e), whose Taylor coefficients at origin is given from the e expansion of ^ as 

B^(t) = 1 - |3C - ^(1 - In 2)j ? + ■ • • . (7.2) 

In order to perform the integration over t in Eq. (17.11) . the analytic continuation of the Borel 
transform B^(t) to the real positive axis of t is necessary. Here we employ the Pade approximant, 
where B^{t) is replaced by the following rational functions 

_ 1+Pi?+--- + Pm^^ 

From Eq. (17.21) . we require that the Pade approximants satisfy pi -qi = -3C+ 1(1 -In 2). Furthermore, 
we incorporate the results around two spatial dimensions in Eq. (|6.32l) by imposing 

^(2-6) = l-6 + --- (7.4) 
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0.5- 




4 



d 



Figure 7.1: The universal parameter ^ as a function of the spatial dimensions d. The upper solid curve 
is from the expansion around d = 4 in Eq. (|3.43l) . while the lower solid line is from the expansion 
around d = 2 in Eq. (16.321) . The middle three curves show the different Borel-Pade approximants 
connecting the two expansions. The diamond indicates the result ^ ~ 0.42 from the Monte Carlo 
simulation f86l|. 



on the Pade approximants as a boundary condition. Since we have three known coefficients from the 
two expansions, the Fade approximants [M/N] satisfying M + N = 3 are possible. Since we could 
not find a solution satisfying the boundary condition ^(2 - e) = I - e for [M/N] = [2/1], we adopt 
other three Fade approximants with [M/N] = [3/0], [1/2], [0/3], whose coefficients pm and qn are 
determined uniquely by the above conditions. 



Fig. 17. II shows the universal parameter ^ as a function of the spatial dimensions d. The middle 
three curves show ^ in the different Fade approximants connecting the two expansions around d = 4 
and d = 2. These Borel-Fade approximants give ^ = 0.391, 0.364, and 0.378 at J = 3, which are 
small compared to the naive extrapolation of the e expansion to J = 3 ^ 0.475). 

The Fade approximant above, however, almost certainly needs serious modification, since it is 
known that in the expansion of ^{e), there exist non-analytic terms at sufficiently higher orders [I102L 



An understanding of the structure of high-order terms in the perturbation theory around d = 4 is 
currently lacking. 




(7.5) 



Chapter 8 

Thermodynamics below Tq 



Now we investigate the thermodynamics of the Fermi gas at finite temperature near the unitarity limit. 
At zero temperature, we found that there exist two difference energy scales in the system; the scale 
of condensate 0o and that of chemical potential yu ~ 60o ^ 0o- Accordingly, we can consider two 
temperature regions where the unitary Fermi gas exhibits different thermodynamics. 

One is the low temperature region where T ~ e0o. In this region, the energy gap of the fermion 
quasiparticle A ~ 0o is still large compared to the temperature. Therefore, thermal excitations of the 
fermion quasiparticle are exponentially suppressed by a factor e'^^'^ ~ e"'^*^. The thermodynamics 
in this region is dominated by the bosonic phonon excitations. The other temperature region is the 
high temperature region where T ~ cpQ. (0o represents the condensate at zero temperature.) In this 
region, the condensate decreases and eventually vanishes at the critical temperature T^. Fermions and 
bosons are equally excited here. We defer our discussion on the high temperature region to Chap. |9] 
and concentrate on the thermodynamics at the low temperature region T «; Tc in this Chapter. 



8.1 Finite temperature formalism 



The extension to finite temperature T follows from the prescription of the imaginary time formalism. 
The system under consideration is described by the following Lagrangian density with the imaginary 
timeO < r < l/T: 



£o = \ dr - - cr-00 1 + (f Wr - — I ^+—^ 



Co 



Co 



Co 



(8.1) 
(8.2) 



(8.3) 



The propagators of fermion and boson are generated by Xo- The fermion propagator is a 2 x 2 matrix. 



G{iCOn,P) 



1 



iOJn + £p -00 
^ -00 i(^n ~ £pj 



(8.4) 
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/iB = ---^--^--^ 

Mb 

diagonal parts : 



off-diagonal part : 




Figure 8.1: Boson's self-energies contributing to the order 0(e). Solid (dotted) lines represent the 
fermion (boson) propagator -G {-D), while the cross in the first diagram represents the //b insertion 
to the boson propagator. The vertex Ho from X2 needs to be added to the second diagram. The last 
diagram gives the off-diagonal part of the self-energy. 



where Sp = p^/2m, Ep = + (p^ and 0o is the condensate in the superfiuid ground state as before. 
The boson propagator D is 

D(/y„,p) = (/y„-^) . (8.5) 

(x>n = 2nT{n + ^) and v„ = InTn are discrete Matsubara frequencies for fermion and boson with an 
integer n = 0, +1, +2, • • • . The unitary Fermi gas around four spatial dimensions is described by the 
weakly-interacting system of fermionic and bosonic quasiparticles, whose coupling g ~ 6^''^ in Xi 



was introduced in Eq. (13.31) as 

/o_2,xl/2 



, ^ , ■ (8.6) 

We define the counter vertex in X2 for the boson propagator by Hoipo, p) = po - Sp/2 and the boson 
chemical potential by /iB = 2// - g^/co. When cq is negative, -g^ Icq - £b gives the binding energy of 
boson to the leading order in e. We consider the vicinity of the unitary point where eb ~ f^o- 

The power counting rule of e developed at zero temperature in Sec. l3.2l holds in the low tempera- 
ture region where the condensate is still large compared to the chemical potential p-jcpo ~ e, while it 
breaks down near the critical temperature because (/iq ^ at T ^ Tc. In the high temperature region 
r ~ Tc, a minor modification of the power counting rule is necessary as we discuss in Chap. HI 



8.2 Phonon spectrum 

The thermodynamics at the low temperature region T «; 00 is dominated by the phonon excitations. 
In order to determine the phonon spectrum, we first study the boson self-energy at zero temperature. 
To the order of 0(e), there are three types of contributions to the boson self-energy as depicted in 



8.2. Phonon spectrum 
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Fig. l8.1[ In addition to the chemical potential insertion //b = + the one-loop diagrams contribute 
to the diagonal part Iln and ofF-diagonal part Tin of the boson self-energy. The vertex Ho from X2 is 
necessary for Iln according to the power counting rule described in Sec. 13. 2[ Then the diagonal part 
of the boson self-energy is given by 



where Ho is defined in Eq. (13.101) and Ha is given by 



(8.7) 



T C dk 



1 



X 



(8.8) 



Ek^p +Ek+p - po 



Ek-E- +Ek+p+ Po 



Since we are interested in physics at the scale of temperature T «; ^q, it is sufficient to evaluate the 
self-energy when the external momentum is small p ~ T (pQ. Expanding Tinip) in terms of pl4>Q 
and performing the k integration with the use of the formula 

f-' T{a)Y{fi - a) 



Jr*co 
*(iTT7 



r(yS) 



we obtain 



nii(p) 



_ 2 r dk El 
^ J {2ny 4 



dk El + sl 3 2^ 
AEl 2 ^ 

k 



(8.9) 



(8.10) 



Similarly, the off-diagonal part of the boson-self energy is given by 



2 fidko 

- V f— 



dk 
dk 



'Pa 



-f) 



y AEk-pEk^p 



1 



-I- 



1 



Ek-Vr + Ek+p - Po Ek-p + Ek+p + Po 



(8.11) 



Expanding Iin{p) in terms of p/^o and performing the integration over k, we obtain 



ni2(p) - g- 



r dk 

J (2^ 



dk (p 



1 

— t 
2 



(8.12) 



As a result of the resummation of these self-energies, the resumed boson propagator D is ex- 
pressed by the following 2x2 matrix: 



i)(po,p) = 



-U21 Di-py'+pB 



n 



22 



(8.13) 



where 1122 = Hn and ni2 = 1121. The dispersion relation of the boson cophip) can be obtained by 
solving the equation det[D'^(a),p)] = in terms of o) as 



(8.14) 
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II 



Figure 8.2: One-loop diagram of boson contributing the effective potential at finite temperature. The 
dotted double line represents the resumed boson propagator T> in Eq. (|8.13l) . 



Note that this expression is valid as long as 0o because of the expansions made to evaluate the 
boson self-energies. Substituting the leading order solution of the gap equation at zero temperature in 
Eq. (13.381) . /zb = ^<Pq, the phonon spectrum is determined to be 

o^Av) = ^y(y + 6'^o). (8.15) 

For the small momentum Sp ecpo, the dispersion relation becomes linear in the momentum as 
'^ph - Cslpl, remaining gapless in accordance with the Nambu-Goldstone theorem. The sound velocity 
of phonon Cs is given by 




(8.16) 



For the large momentum Sp » 60o: the dispersion relation approaches that of the free boson as 

<^ph - £p/2. 



8.3 Effective potential and condensate 

At finite temperature, the phonon excitations contribute to the effective potential, and consequently, 
the magnitude of the condensate decreases. The temperature dependent part of the effective potential 
Vt{(Pq) to the lowest order in 6 is given by the one-loop diagram of the boson with the resumed 
propagator in Eq. (|8.13l) [Fig. 18.211 : 



Wo) = |j] r-^Trln[D(/v„,p)-i] 

" (8.17) 

(2ny L J 



The n-th interaction vertex among phonons (p" is of the order e"^^ and appears in the effective potential 
only at higher orders. Then the contribution of VT((f>o) to the gap equation is 

— = 7T^^B(^^ph)^— , (8.18) 



8.4. Thermodynamic functions at low temperature 
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where fB{x) = 1 /{e^^^ - 1) is the Bose distribution function and da)ph/d(f)o is given from Eq. (18.141) by 

dconh 3eo + 2600 
d(f>o 

There are two limiting cases where the integration over p in Eq. (|8.18l) can be analytically per- 
formed. Since the integral is dominated by the integration region where Sp ~ T, we can approximate 
the phonon spectrum by its linear branch ojphip) - Cs\p\ when the temperature is very low T «; ecpo- 
In this case, the integration over p in Eq. (18.181) leads to 

SVM^mLl:^.)', (8.20) 

d(pQ (po \2n I 

On the other hand, when the temperature is located in the intermediate region e0o ^ T (p^, the 
phonon spectrum can be approximated by its quadratic branch oj^^ - Spjl. In this case, the integration 
over p in Eq. (|8.18l) results in 

d(po 4 

Now, from the gap equation 5(Veff(0) + yT{(p))/d(p = with (p = (f>o + (f>T, one finds the temperature 
dependent correction of the condensate (pj satisfies 

where Veff is the eff'ective potential at zero temperature in Eq. (13.361) . To the leading order in e, (pr at 
T «; 600 is given by 

(pT = — , (8.23) 

while at 60o T ^ (pQ, 

(pT = -e-—. (8.24) 

(po 

The condensate in total is = 0o + 0r, which decreases as the temperature increases. Note that since 
cpT 600, the leading part of the condensate does not change in the temperature region considered 
here T «; 0o. The effective potential is given by the sum of the zero temperature and finite temperature 
parts; VeffC^o + (pr) + W(0o + 0r) - V^eff(0o) + Vt{(Po)- 



8.4 Thermodynamic functions at low temperature 

The temperature dependent part of the pressure Pph at the low temperature region T «; 0o is given 
from the eff'ective potential in Eq. (|8.17l) by 

/'ph = - Vr(0o) = - J-^T\n[l- ^,-p^(P)/^] . (8.25) 
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The phonon contributions to the fermion number density, the entropy density, and the energy density 
are computed from the thermodynamic relations, A^ph = dPph/d/u, 5ph = dPph/dT, and Eph = {uNph + 
rSph - Pph, respectively. Here we show analytic expressions for these thermodynamic functions in 
the two cases where the analytic evaluation of the p integration in Eq. (|8.25l) is available. 

When the temperature is very low T «; ecpo, only the linear branch of the phonon spectrum 
^ph(p) - Cslpl is important to the thermodynamic functions. In this case, the integration over p can 
be performed analytically to lead to 



■^Ph - 4 ~ 2 — r J, \2- (O./D) 



Accordingly, we obtain the phonon contributions to the fermion number density, the entropy density, 
and the energy density; 



48^(5) m^T^ 



A^ph = ^TTTTT' (8-27) 



(e(/)oy 

24^(5) m'T' /, ^ fib \ 
Epb = — i— T-^a l + -r • (8-29) 
(60o)^ \ e0o/ 

Since actual experiments or simulations are performed with the fixed fermion density, it is useful 
to show the thermodynamic functions at fixed instead of fixed /i. From Eqs. (13.381) . (|3.40l) . and 
(18.271) . we find the chemical potential for the fixed fermion density increases as a function of the 
temperature as 

//= ;Uo + 48^(5)— , (8.30) 

where yUo represents the chemical potential at zero temperature in Eq. (|3.42|) . Normalizing // by the 
Fermi energy in Eq. (I3.41|) . we have 

E^ti.mim', (8.31) 

The other thermodynamic functions for the fixed fermion number density are given by 

1 , (8.32) 

(8.33) 
(8.34) 



p 




3^(5) 
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8.4. Thermodynamic functions at low temperature 
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where Pq and Eq represent the pressure and energy density at zero temperature in Eqs. (|3.44l) and 
(13.451) . respectively. These expressions are valid in the low temperature region where T 60o- 

On the other hand, when the temperature is located in the intermediate region ecpo «; T «; (pQ, we 
can expand the phonon spectrum tOphip) in terms of ecpo/sp up to its first order 

^ph(p) - y + (8-35) 
In this case, the integration over p can be performed analytically again to result in 

^(3) , , m^T^ 
/'ph - ^m^T' - -p^60o. (8.36) 
12 

Accordingly, we obtain the temperature dependent parts of the fermion number density, the entropy 
density, and the energy density; 



6 ' 



A^ph = —, (8.37) 



3^(3) 2rr2 rn^T 

Sph = — 5— m T —e(po, (8.38) 

6 



2/(3) , , m^T^ I £b \ 

£ph = -^m^T' - —^ecp, 1 - . (8.39) 

From Eqs. (13.381) . (|3.40l) . and (|8.37l) . we find the chemical potential for the fixed fermion density 
increases as a function of the temperature as 

+ 6^3-- (8.40) 
6 00 

Normalizing // by the Fermi energy in Eq. (13.411) . we have 

/i = ^+,3/2^(n'. (8.41) 



Sp £f 6 \£f 

The other thermodynamic functions for the fixed fermion number density are given by 



eyiN spN \£p) 6 yep 



E Eq ^ ^ (TV .ijTi'^lT 
-^ + 8Z-(3) — -e^'^ 



SfiN SfN yep/ 3 yep 



2 



+ 8^(3) - - ' (8-43) 



^ = 12^(3)1-1 -6^/'=-— . (8.44) 
N \ep/ 3 ep 

These expressions are valid in the intermediate temperature region where €(po «; T «; 0o- 
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— 



0.5- 




Figure 8.3: The behavior of the condensate as a function of the temperature T to the lowest order 
in 6. (po = jUB/e is the value of the condensate at T = and the critical temperature is located at 
Tc/(f>o = 1/(2 In 2) « 0.721348. 



8.5 Effective potential near Tc 



At the end of this Chapter, we study the behavior of the condensate as a function of the temperature 
for a given fi. The effective potential to the leading order in e is given by one-loop diagrams of 
fermion with and without one fi insertion, which is equivalent to the mean-field approximation. Since 
the critical dimension of the superfluid-normal phase transition is four, the mean-field approximation 
remains as a leading part at any temperature in the limit d ^ 4. Then the leading contribution to the 
effective potential at finite temperature is given by 

dp 



Ep-^fi + 2rin(l + e-^-'^) + MEp)^fi 



where f^ix) = 1 /(e'^^^ -i- 1) is the Fermi distribution function. 

For the low temperature T «c 0, we can neglect the exponentially small factor e'^p^^ «; 1 and the 
integration over p reproduces the leading order effective potential at zero temperature in Eq. (|3.36l) . 
In the opposite limit where T, we can expand Veff in terms of (p/T to lead to 



Veff(0) = yeff(O) + 



rin2-^ 
2e 



In) \6T\2nl 



(8.45) 



where /zb = 2fi + Sb. From the coefficient of the quadratic term in 0, we can read the critical tempera- 
ture Tc to the leading order in e as 



T = 



e2\n2 

and the value of the condensate just below T^ as 



+ 0(6), 



= 8r(rc -r)in2 + c>(6). 



(8.46) 



(8.47) 



The critical exponent of the condensate (p ~ {T^ - T)^^^ will be shifted if we include higher order 
corrections to the effective potential. The condensate (p in the intermediate range of the temperature 
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is obtained by solving the gap equation dVes/d^) = 0; 

dSp^MEp) = 0. (8.48) 
The numerical solution of the gap equation as a function of T is shown in Fig. 18.31 
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Chapter 9 

Thermodynamics above Tq 



9.1 Power counting rule of e near 

In the e expansion at the zero or low temperature region, the chemical potential is small compared to 
the condensate // ~ 60o and we made expansions in terms of /i/0o ~ f as well as the small coupling 
g ~ e^l'^. Near the critical temperature, the ratio ///(/> is no longer small because (p vanishes at T = Tc, 
but jJi/Tc is 0{e) as it is clear from = fx^/iellnl). Therefore, we can still treat the chemical 
potential as a small perturbation near and the same power counting rule of e described in Sec. 13.21 
holds even above Tc just by replacing with T. Hereafter we consider T ~ T^to be 0(1). 



9.2 Boson's thermal mass 

First we study the self-energy of boson at T > Tc. The leading contribution to the self-energy is the 
chemical potential insertion /j.^ as well as the one-loop diagram rin shown in Fig. l8.1[ 



(9.1) 



2 r dk 

2 r dk I - fA^k-p/i) - fviek+pii) 
" J (2^ 2sk - iv + ep/2 " 

For the zero Matsubara frequency mode v„ = at the small momentum Sp <^ T,we have 

nn(0,0)-/ r^^ = ,r21n2. (9.2) 

Therefore, the zero Matsubara frequency mode has the non-negative thermal mass Ut = 6 T 2 In 2 - 
//b ~ 6 at r > Tc. The condition of the vanishing thermal mass Ut = gives the critical temperature 
Tc = /^B/(f 21n2) equivalent to Eq. (18.461) . As we will see below, at a sufficiently high order in the 
perturbation theory near Tc [e^ (e) compared to the leading term in the pressure (fermion density)], 
the resummation of the boson self-energy is needed to avoid infrared singularities appearing in the 
zero Matsubara frequency mode. 
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0(1) 0(1) 0(e) 

Figure 9.1: Three types of diagrams contributing to the pressure up to the next-to-leading order in e. 
Each /i (//b) insertion to the fermion (boson) line reduces the power of e by one. 



9.3 Pressure 



Now we calculate the thermodynamic functions at T > to the leading and next-to-leading orders in 
e. There are three types of diagrams contributing to the pressure up to the next-to-leading order in e as 
depicted in Fig. 19.11 one-loop diagrams with and without one /i {p.^) insertion and two-loop diagram 
with a boson exchange. Note that at T > T^, the boson's one-loop diagram contributes as 0(1) as well 
as the fermion's one-loop diagram. Then the pressure from the one-loop diagrams is given by 



Pi= [2r In (l + e-'-l^) - r In (l - e-'-'^^) + l/uMsp) + 



= T 



11 9 In 2 ^(3) +11^(3) ^n^M ^2n^iuj, 
y^(3) ^ e + + — y 



The contribution from the two-loop diagram to the pressure, which is 0(e), is given by 
Pi 



2 2 V Cdpdq 

> j j:^--^Gn{mi,p)G22(ioj,„,q)Diici)n- ioJni,p- q) 



n,m 



/dpi 
(2^ 



yd 



Ep "1- Sq Sp~ql'2. 



The numerical integrations over p and q result in 



(9.3) 



(9.4) 



(9.5) 



where Cp ~ 8.4144. From Eqs. (19.31) and (19.51) . we obtain the pressure up to the next-to-leading order 
in e as 



P = P,+P2 

11 



= T 



91n2^(3)+lir(3) 
^^(3) ^ e 



eCp + —— + 



6 T 3 T 



(9.6) 



The entropy density S and the energy density E to the same order can be computed from the thermo- 
dynamic relations S = dP/dT and E = f^N + TS - P. 
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9.4 Fermion number density 



The fermion number density to the next-to-leading order in e can not be obtained simply by differenti- 
ating the pressure in Eq. (19.61 ) with respect to the chemical potential A'^ = dPjdix. Since the pressure to 
the leading order in e does not depend on jj. and the yu derivative djdiJ. ~ 1 / e enhances the power of e 
by one, we need to compute the one-loop diagrams with two ju (jUb) insertions and two-loop diagrams 
with one ip.^) insertion. Then the fermion density from the fermion's one-loop diagrams is given by 



J w 



dp 

2 „1 



nnn2 + 6C(.2) 2 In 2 
6 -2 '^—^ 



lmT\ 
I 2^/ 



(9.7) 



On the other hand, the boson's one-loop diagrams contribute to the fermion density as 



=2 



dp 

Y 



Uspl2) - ^/B(£p/2)/B(-ep/2) 



(9.8) 



Apparently, the last term has an infrared singularity because the Bose distribution function behaves 
as fB{£p/2) - 2T/Sp at the small momentum Sp T. In order to resolve this infrared singularity, 
the resummation of the boson self-energy Hn = 21n2 is needed at the small momentum region 
Sp ~ ju. Within the accuracy we are working, we can rewrite A^b as 

A^B = 2 J ^ [/B(ep/2 + n^) - 6 2 In 2 Msp/2)M-ep/2)] , (9.9) 

where Tlj = eT 2\n2 - jub. Now the first term is infrared finite, where the boson's thermal mass 
^t/T = (2 \n2)e- ubIT plays a role of an infrared cutoff". Integrating over p and expanding up to the 
next-to-leading order in e, we have 



An'^ 27r2ln2+ 12^'(2) 
~i 3 



6-8(l-ln^l^ + 0(62) 



I 2^/ 



(9.10) 



The logarithmic term ~ InTlr/T appears as a consequence of the resummation. The second term 
in Eq. (|9.9I) . which is still infrared divergent, will cancel with the infrared singularity existing in the 
two-loop diagram. 

The contribution from the two-loop diagrams to the fermion density is given by 



A^9 



, Cdkt 



dp 

"yd 



(9.11) 



X 



/F(gfc+p/2)/F(-gfc+p/2) [/F(gfc-p/2) + MSp/2)] - 2/F(gfc+p/2)/B(gp/2)/B(-gp/2) 
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The second term in the numerator contains the infrared singularity at small Sp. Extracting the diver- 
gent part, we can rewrite A^2 as 

^dkdp /F(efc+p/2)/F(-efc+p/2) [fFi^k-p/l) + /B(ep/2)] 

2 I ^ dk dp 2 [fvisk+p/i) - /F(gfc)] /B(gp/2)/B(-gp/2) 



N, 



2 Cdka 
J (2^ 

/dka 

/dkt 



dp fvisk) 



fB(Sp/2)M-ep/2). 



(9.12) 



rr' SkT 

One finds the k integration in the last term can be performed to lead to Hn = 6 T 2 In 2 in Eq. (19.21) . 
which exactly cancels out the infrared divergent part in Eq. (19.91) . The numerical integrations over k 
and p in the first two terms result in 



A^2 = -CNe 



\2nl " ""'"J {2ny 



641n2 



J (2;r 



fBiSp/2)M-Sp/2), 



(9.13) 



where Cn ~ 1.92181. Gathering up all contributions, Eqs. (19.71) . (19.91) . and (19.131) . the fermion number 
density to the leading and next-to-leading orders is given by 



Sn^ 3;r2ln2 + 18^'(2) 



6- 6Cm + 



2 4 ^ 

We define the Fermi energy ep through the relationship in Eq. (13.411) as 



\2^) 



(9.14) 



2n 


1 






m 








2 



2 
1 + 

-I- 



27-3-21n2 3^'(2) Cn 1, 

+ - In 



8 

21n2/z 8 
In^T ~ 3^ 



2n^ 8 

Zi 



37r^ 
2 



(9.15) 



The logarithmic correction (Jij IT) In Ilr/r ~ e In 6 is a consequence of the resummation to avoid the 
infrared singularities, while it vanishes at the critical temperature IIt-^ = 6 Tc 2 In 2 - /^b = 0. 



9.5 Critical temperature 



The critical temperature in units of the Fermi energy directly follows from Eq. (19.151) with the use of 
the relationship 2/z -I- Sb = f T'c 2 In 2: 



Sp 



2 

3^ 



,2y-3-21n2 3^'(2) Cw-2(ln2)2 1 (3n^ 

'-^ s s^^^sH—^'' 



+ 



ln2 eb 
3n^ 



= 0.260 - 0.01 12 6 + 0.0234 — + 0(e\ 



(9.16) 
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where the numerical value Cn ~ 1.92181 is substituted. We find the critical temperature Tc is an 
increasing function of the binding energy eb near the unitarity limit. The next-to-leading order cor- 
rection is reasonably small compared to the leading term even at 6 = 1 . The naive extrapolation of 
the critical temperature to the physical case of J = 3 gives Tc/sf ~ 0.249 in the unitarity limit Sb = 0. 
This value is surprisingly close to results from two Monte Carlo simulations, Tdsp = 0.23(2) lf89l 
and Tc/sp « 0.25 [|92l . while other two simulations provide smaller values, Tdsp < 0.14 [|90l and 
rc/eF = 0.152(7)11911. 

It is also interesting to compare the critical temperature in the unitarity limit with that in the BEC 
limit Tbec- In the BEC limit, all fermion pairs are confined into tightly bound molecules and the 
system becomes a non-interacting Bose gas where the boson mass is 2m and the boson density is 
N/2. The critical temperature for the Bose-Einstein condensation of such an ideal Bose gas at J > 2 
spatial dimensions becomes 

T^BEC 1 



2 2 



2/d 

(9.17) 



To the leading and next-to-leading orders in e = 4 - d, the ratio of the critical temperatures in the 
unitarity limit Tc and in the BEC limit Tbec at the same fermion density is given by 



BEC 



- [l + 0.0177 6 + 0{e^)] = 0.943 + 0.0167 e + 0{e^). (9.18) 



The ratio is slightly below unity, indicating the lower critical temperature in the unitarity limit Tc < 
Tbec- The leading order term of the above ratio, Tc/T^ec = V8/9, has the following clear physical 
interpretation: The critical temperature for the Bose-Einstein condensation at J = 4 is proportional 
to a square root of the boson's density. In the BEC limit, all fermion pairs form the bound bosons, 
while only 8 of 9 fermion pairs form the bosons and 1 of 9 fermion pairs is dissociated in the unitarity 
limit [see the leading order terms in Eqs. (19.71 ) and (19.101 )1. Thus, their ratio in the critical temperature 
should be Tc/Tqec = V879 < 1 at J = 4. 

The more appropriate estimate of r^. at J = 3 will be obtained by matching the e expansion with 
the exact result around d = 2. The critical temperature at unitarity in the expansion over e = d - 2 
is given by Tc = (e^/^r) A, where A/sp = (2/e) e'^^'^ is the energy gap of the fermion quasiparticle at 
zero temperature [|94]| : 

T 2e'^~^ 

e-^l'\_l + 0{e)\. (9.19) 



We shall write the power series of e in the form of the Borel transformation. 



Tc{e) 2ey-' 



ep 



— e-^" dte-'Briet). (9.20) 
^ Jo 



Bj^^t) is the Borel transform of the power series in Tde), whose Taylor coefficients at origin is given 
by firJO = 1 -I- • • • . In order to perform the integration over t in Eq. (|9.20l) . the analytic continuation 
of the Borel transform firJO to the real positive axis of t is necessary. Here we employ the Pade 
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Tcle^ 
0.3| ' 




Figure 9.2: The critical temperature Tc at unitarity as a function of the spatial dimensions d. The 
upper solid line is from the expansion around d = Am Eq. (19.161 ). while the lower solid line is from 
the expansion around J = 2 in Eq. (19.191 ). The middle two curves show the different Borel-Pade 
approximants connecting the two expansions. The symbols at J = 3 indicate the results from the 
Monte Carlo simulations; TJs^ = 0.23(2) [[891 (circle), TJsf < O.U ^ (down arrow), TJsf = 
0.152(7) mi (diamond), and TJsf « 0.25 |l92l (square). 



approximant, where BjXt) is replaced by the following rational functions 

BrM}= ' (9.21) 

Then we incorporate the results around four spatial dimensions in Eq. (19.161) by imposing 

^^^^^^^ =0.260-0.01 12e+--- (9.22) 

on the Pade approximants as a boundary condition. Since we have two known coefficients from the 
6 expansion, the Fade approximants [M/A^] satisfying M + = 2 are possible. Since we could not 
find a solution satisfying the boundary condition in Eq. (19.221 ) for \MIN^ = [1/1], we adopt other two 
Fade approximants with [M/N] = [2/0], [0/2], whose coefficients p,„ and qn are determined uniquely 
by the above conditions. 

Fig. l9.2l shows the critical temperature in units of the Fermi energy as a function of the spatial 
dimensions d. The middle two curves show Tc/sf in the different Fade approximants connecting the 
two expansions around d = 4 and d = 2. These Borel-Fade approximants give Tc/s^ = 0.173 and 
0.192 at J = 3, which are located between the naive extrapolations to J = 3 from the e = 4- - d 
expansion (Tc/sf 0.249) and the e = d - 2 expansion {Tc/sf 0.153). It is also interesting to 
compare our results with those from the recent Monte Carlo simulations, where Tc/sf = 0.23(2) [89], 
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TJsf < 0.14 [I90l, TJsf = 0.152(7) and TJsf ~ 0.25 Although these resuks from the 
Monte Carlo simulations seem not to be settled, the interpolation of the two expansions provides the 
moderate value Tc/ep = 0.183 + 0.014 not too far from the Monte Carlo simulations. 



9.6 Thermodynamic functions at Tc 

Finally we show the thermodynamic functions at in the unitarity limit £b = to the leading and 
next-to-leading orders in e. The pressure P normalized by the fermion density s^iN follows from 
Eqs. (|9^ . (19.141) . and (19.151) . Introducing the numerical values Cp ~ 8.4144 and Cn ~ 1.92181, we 
obtain the pressure up to the next-to-leading order in e as 



= 0.116 + 0.0188 6. 



(9.23) 



From the universal relationship in the unitarity limit E = (d/2)P, the energy density is given by 

E 



= 0.232 - 0.0205 e. 



(9.24) 



The chemical potential at the critical temperature fi = 6 Tc In 2 is 0(e). Normalizing /i by the Fermi 
energy in Eq. (|9.16l) . we have 



= eln2 



3n^ 



= 0.180 6. 



Then the entropy density TcS = (d/2 + 1)P - /uN is given by 

= 1.340- 0.642 e. 



S_ 

N 



(9.25) 



(9.26) 



The next-to-leading order corrections to the pressure and energy density are reasonably small com- 
pared to the leading order terms, while that is large for the entropy density. 

We match the thermodynamic functions at in the expansions over e = 4 - d with those around 
J = 2 as we demonstrated for T^/sf. The critical temperature around d = lis Tc/sf ~ e"'''*^, which is 
exponentially small and negligible compared to any power series of e. Therefore, the pressure, energy 
density, and chemical potential at T^. in the expansions over 6 = J - 2 is simply given by those at zero 
temperature [94 1; 



E 



E. 



Tc 



2 1 5_ 

d + 2^ 2 8^' 



d 1 3_ 
d + 2^ ~ 2~ 8^' 



^ = 1 - 6- 



(9.27) 
(9.28) 
(9.29) 
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Tc/sf 




E/(£fN) 




S/N 


6 expansion (NLO) 


0.249 


0.135 


0.212 


0.180 


0.698 


Borel-Pade approximant 


0.183 


0.172 


0.270 


0.294 


0.642 


self-consistent approach II78II 


0.160 


0.204 


0.304 


0.394 


0.71 


Monte Carlo simulation If89]| 


0.23(2) 


0.27 


0.41 


0.45 


0.99 


Monte Carlo simulation 11911 


0.152(7) 


0.207(7) 


0.31(1) 


0.493(14) 


0.16(2) 



Table 9. 1 : Comparison of the results by the e expansion with other analytic and numerical calculations 
in the unitarity limit. 



4Zl = ^,-./^ ,9.30) 



A straightforward calculation shows that the entropy per particle at to the leading order in e is 
given by 

5 

Using the Borel-Pade approximants connecting the two expansions above, thermodynamic func- 
tions at J = 3 are found to be P/isFN)\T^ = 0.172 + 0.022, E/ispN)\T^ = 0.270 + 0.004, ///splr, = 
0.294 + 0.013, and S /N\y^ = 0.642. The errors here indicates only the uncertainty due to the choice 
of different Fade approximants. In Table 19.1 1 we summarize our results on the critical temperature 
and the thermodynamic functions at Tc as well as results from other analytic and numerical calcula- 
tions. In the recent Monte Carlo simulation, the thermodynamic functions at the critical temperature 
is P/(£fN)\t^ = 0.207(7), E/(£fN)\t^ = 0.31(1), yu/splr, = 0.493(14), and S /N\t^ = 0.16(2) M^- We 
see that the interpolations of the two expansions indeed improve the series summations compared to 
the naive extrapolations from d = 4 or d = 2, while there still exist deviations between our results and 
the Monte Carlo simulation. We can understand these deviations partially due to the difference in the 
determined critical temperature. The large deviations existing in ///sp and S/N may be because we 
know only the leading term for /i and the next-to-leading order correction for S is sizable. 



Chapter 10 

Summary and concluding remarks 



We have developed the systematic expansion for the Fermi gas near the unitarity limit treating the 
dimensionality of space d as close to four. To the leading and next-to-leading orders in the expansion 
over 6 = 4 - d, the thermodynamic functions and the fermion quasiparticle spectrum were calculated 
as functions of the binding energy S], of the two-body state. Results for the physical case of three 
spatial dimensions were obtained by extrapolating the series expansions to e = 1. We found the 
universal parameter of the unitary Fermi gas to be 

^= — = \e^'^ + ^e^'^ In 6 - 0.0246 e^'^ + 0{e^'^) « 0.475. (10.1) 
£f 2 16 

The fermion quasiparticle spectrum in the unitarity limit was given by the form a»F(p) - yj{£p - sqY + 
with the energy gap 

- = - -0.691 + 0(6) « 1.31 (10.2) 

J" 6 

and the location of the minimum of the dispersion curve 

— = 2 + 0(6) ~ 2. (10.3) 
J" 

We also found the condensate fraction in the fermion density to be 

^ = 1 - 0.0966 6 - 0.2423 + 0{^) « 0.661. (10.4) 

Although we have only the first two terms in the e expansion, these extrapolated values give reason- 
able results consistent with the Monte Carlo simulations or the experimental measurements. Further- 
more, the corrections are not too large even when extrapolated to 6 = 1, which suggests that the picture 
of the unitary Fermi gas as a collection of weakly interacting fermionic and bosonic quasiparticles 
may be a useful starting point even in three spatial dimensions. 

We have also formulated the systematic expansion for the unitary Fermi gas around two spatial 
dimensions. We used the results around J = 2 as boundary conditions which should be satisfied 
by the series summations of the expansion over e = A - d. The simple Borel-Pade approximants 
connecting the two expansions yielded ^ = 0.378 + 0.013 at <i = 3, which is small compared to 
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the naive extrapolation of the e expansion [Fig. 17.111 . In order for the accurate determination of ^ 
at d = 3, the precise knowledge on the large order behavior of the expansion around four spatial 
dimensions as well as the calculation of higher order corrections would be desirable. Once these 
information become available, a conformal mapping technique, if applicable, will further improve the 
series summations [1133L 

In Chapters |4] and |5l the phase structure of the polarized Fermi gas in the unitary regime with equal 
and unequal masses has been studied based on the e expansion. The gapless superfluid phase and the 
superfluid phase with spatially varying condensate were found to exist between the gapped superfluid 
phase and the polarized normal phase in a certain range of the binding energy and the mass difference. 
In the equal mass limit, our study gives a microscopic foundation to the phase structure around the 
splitting point which has been proposed using the effective field theory [|112| . Moreover, we found the 
gapless phase with spatially varying condensate around d = 4 exists only at 0.494 < (akp)'^ < 1 and 
terminates near the unitarity limit. Our result suggests that the phase with spatially varying condensate 
existing in the unitary regime may be separated from the FFLO phase in the BCS regime ni24lll25ll . 
We also found the splitting point and the same phase structure around it in the unitarity limit but with 
finite mass difference when the majority is the lighter fermions. The range of the superfluid phase 
with spatially varying condensate can be estimated to be 1.37 < mminority/''^majority ^ 3. The splitting 
points form a smooth line in the three-dimensional phase diagram of 1 l{akp), the mass difference 
K, and polarization chemical potential H. It gets away from the unitarity with increasing the mass 
of major fermions. Eventually the splitting point becomes unstable due to the competition with the 
polarized normal state at the point given by {ak^Y^ ~ 1.81 and mmajority/'^mmority * 5.54. Further 
investigation will be worthwhile to confirm these possibility. 

In Chapters [8] and [9l the thermodynamics of the Fermi gas near the unitarity limit at finite tem- 
perature has been investigated using the systematic expansion over e = A - d. We discussed that the 
thermodynamics in the low temperature region T Tc is dominated by the bosonic phonon excita- 
tions. The analytic formulas for the thermodynamic functions at the fixed fermion density are derived 
in the two limiting cases; T «; eT^ in Eqs. (IOT]) - (l834l) and eT^ ^ T ^ in Eqs. (l84T]) - (l8l4l) . 
In the high temperature region T ~ T^, the fermionic quasiparticles are excited as well as the bosonic 
quasiparticles. We showed that the similar power counting rule of 6 to that developed at zero temper- 
ature works even above T^. The critical temperature and the thermodynamic functions around T^. 
were calculated to the leading and next-to-leading orders in e. We found the critical temperature is an 
increasing function of the binding energy near the unitarity limit: 



The next-to-leading order correction is reasonably small compared to the leading term even at e = 1 . 
In the unitarity limit = 0, the naive extrapolation of the critical temperature to the physical case of 
d = 3 gives TJe^ ~ 0.249. 

We also discussed the matching of the e expansion with the expansion around d = 2. The critical 
temperature at unitarity in the expansion over e = J - 2 is given by Tc/sf = (2e^~^ /n) e'^^'^, whose 



— = 0.260 - 0.01 12 6 + 0.0234 — + O(e^). 



(10.5) 
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Figure 10. 1 : The Fermi gas at infinity scattering length akp ^ oo as a function of the spatial dimension 
d. At d < 2, such a system becomes a free Fermi gas, while the limit of zero-range interaction 
rokp can not be taken at J > 4. Therefore, the non-trivial universal limit exists only at 2 < J < 4. 



naive extrapolation to 6 = 1 gives Tc/ep * 0.153. The Borel-Pade approximants connecting the 
two expansions yielded Tc/ep = 0.183 + 0.014 at J = 3, which is a moderate value located between 
the two naive extrapolations [Fig. 19.211 . These values are not too far from the results obtained by 
the recent Monte Carlo simulations where TJsp = 0.15 ~ 0.25 [[891 |9ll |92l|. We also applied the 
Borel-Pade approximants to the thermodynamic functions at T^., which yielded P/(sy^N)\j^ ~ 0.172, 
E/(spN)\r^ « 0.270, m/sf\t, « 0.294, and S /N\t^ « 0.642 atd = 3. 

We conclude that the e expansion is not only theoretically interesting but also provides us an 
useful analytical tool to investigate the properties of the unitary Fermi gas. Although our results 
are consistent with those from Monte Carlo simulation, there are several problems which should 
be clarified in future. For example, the Borel-Pade approximants employed here to match the two 
expansions around four and two spatial dimensions do not take into account the large order behavior 
of the expansions around d = 4 and d = 2. Most probably, the expansions over e = 4-d and e = d-2 
are not convergent because d = 4 and d = 2 corresponds to the quantum phase transition points as 
depicted in Fig. 110.11 Thus the Borel transform B(t) of such series expansions will have singularities 
somewhere in the complex ?-plane. Furthermore, an understanding of the analytic structure of high- 
order terms in the perturbation theory around J = 4 is currently lacking. In order for the accurate 
determination of the universal quantities, e.g., ^, Tdsp, at J = 3, it will be important to appropriately 
take into account the knowledge on the large order behaviors of the expansions both around four and 
two spatial dimensions. The calculation of higher order corrections to our results is also important 
for this purpose. In Appendix, we show a type of diagrams which grows as a factorial of n at order 
e" and may be related with the large order behavior of the e expansion. These problems should be 
further studied in future works. Also the application of the e expansion to investigate other problems 
of the unitary Fermi gas, including its dynamical properties (dynamical response functions and kinetic 
coefficients) and the structure of superfluid vortices, will be extremely interesting. 
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Appendix A 

Large orders in the e = 4 - d expansion 



In this Appendix, we show that there exists a type of diagrams which grows as n! by itself at order e" 
of the e = 4 - d expansion. Such a factorial contribution originates from the large momentum region 
of the loop integrals which resembles the ultraviolet renormalon in relativistic field theories I1137[ 
I138[|139l . An example of the n + 1-loop diagram contributing to the eff"ective potential as n! at 0(6") 
is depicted in Fig. I A. 11 which can be written as 

dk 



^i_ r dk 

" nj (2;r)^+i 



[iUoik) + U,ik))Dik)r. 



(A.l) 



where Tlo(k), Tla(k), and D(k) are respectively defined in Eqs. (13.101) . (13.171) . and (13.81) . Introducing 
these definitions and integrating over k^, we obtain the following expression for V„, 
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Then we consider the q integration in the bracket. Since the q integration contains a logarithmic 
divergence at J = 4, we subtract and add its divergent piece as 
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Appendix A. Large orders in the e = A - d expansion 




Figure A.l: A n-ih. order diagram dXd = A which contributes to the effective potential as n! by itself. 
The counter vertex -illo for each bubble diagram is understood implicitly. 



The q integration in the second line becomes finite at d = A and does not produce singular loga- 
rithmic terms. So we concentrate on the q integration in the last line, which can be evaluated in the 
dimensional regularization as 
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Then V„ in Eq. (IA.2I) to the order t" becomes 
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By picking up the logarithmic terms, we find the p integration from its large momentum region 
Sp » Sq and Sp » 0o gives the following contribution to the efi"ective potential, 

,3 



-(f)V(n), 



(A.6) 
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which grows as a factorial of n. Note that this nl contribution comes from the single diagram depicted 
in Fig. lA.ll At the moment, it is not clear whether such a nl contribution survives at 0(e") of the 
effective potential and dominates the large order behavior of the e expansion. Possibly there might be 
cancellations with other n-th order diagrams and/or with subleading contributions from lower-order 
diagrams. 
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